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Preface 


The inexpensive Timex/ Sinclair 1000 computer is often purchased as a first entry 
into the computer field. But its toy-like appearance discourages many serious pro- 
spective buyers who do not realize its considerable potential. They try a few dem- 
onstration programs, find the computer awkward to use, and put it aside. This 
book, written specifically for the Timex/Sinclair 1000, opens up new fields of 
applications for the high school or college student, and the engineer or scientist. 

A high school or college student will recognize many problems treated here, 
such as complex arithmetic, combinatorics, distributions, and solutions to tran- 
scendental and differential equations, as the subjects of his or her homework as- 
signments. And, the professional will find that with these programs many mathe- 
matical and engineering problems can be easily solved on the little computer. 

Science and Engineering Programs for the Timex/ Sinclair 1000 contains a col- 
lection of progams I developed over a number of years for use in my profession as 
a communications engineer and for use in my hobbies of home computing and 
electronic design. They cover assorted problems in the fields of electrical engineer- 
ing, probability, statistics, queuing theory, reliability, curve fitting, graph genera- 
tion, number theory, computer science, artificial intelligence, and other 
disciplines. 

Some of these programs were originally developed for the Hewlett-Packard 
HP-25 and HP-67 programmable calculators and were later translated into 
BASIC. In the process of translation I added a number of features relating to bet- 
ter prompting and input error trapping. Many more useful programs were written 
in BASIC originally and added to the collection. A few of the programs have been 
published in magazines. Where this is the case a proper acknowledgement to the 
magazine is given. 

The coding here has been developed from earlier, less sophisticated versions de- 
signed to run on the Radio Shack/Sharp Pocket Computers. Translation into 
Timex/Sinclair 1000 BASIC required many changes, as the Sharp and Timex/ 
Sinclair BASIC languages are only partially compatible. Many programs were 
completely rewritten and enhanced to make use of advanced Timex/ Sinclair 1000 
features such as string manipulation and screen display. In the course of writing 
the Timex/Sinclair 1000 programs I developed many practical programming 
hints which are detailed in the Appendix of this book. These hints, and various 
programming techniques described in the Program Remarks section of each 
chapter, should further sharpen the reader's programming skills. 

All programs in Science and Engineering Programs for the Timex/ Sinclair 
1000 are specifically written for the computer with at least 2K of RAM or a 
similarly configured ZX-81. However, each listing is preceded by a program 
description and the underlying equations, which should make translation to other 
BASIC computers quite easy. When translating into other computer languages be 
aware of the idiosyncrasies of each language. For example, the integer function 
INT( — 3.5) will return — 4.0 in BASIC, but it will return — 3.0 in Fortran! 

I hope you will find these programs as useful as I do, and that they will let you 
enjoy your Timex/Sinclair 1000 even more. 

Cass Lewart, Holmdel, NJ, 1983 


Introduction 


Each program chapter consists of five sections: 

1. Program Description. The Program Description provides the reader with the 
necessary background and understanding of the problem. 

2. Instructions. The Instructions demonstrate the key sequences required to 
execute the program on the Timex/ Sinclair 1000. 

3. Examples. Examples further clarify the use of the program. This section in- 
cludes paragraphs on problem definition, discussion of results, and a sample 
listing of the program run directly from the computer by means of the COPY 
command. 

4. Programming Remarks. The remarks explain program flow and program- 
ming techniques. They also substitute for the REMark statements in the listings 
omitted because of tight memory limitations on the 2K computer. 

5. Program Listings. Printed with the LLIST command, and interfaced to a 
standard dot matrix printer, each listing corresponds exactly to the program 
which generated the examples. 

You will notice that all programs contained in this book have certain common 
features. They are easy to use and provide comprehensive prompting and error 
trapping. The programs are not trivial; e.g., they do not do something which 
could be done more easily with a hand calculator. With very few exceptions, each 
program effortlessly performs functions which would be difficult to perform in 
any other way. None of the programs is simply a formula translation. And, they 
do not use undocumented tricks to save a byte of memory here and there at the 
cost of clarity. 


Before keying in the programs you may want to read the Appendix section 
describing various programming hints. And finally, last but not least, read the 
User’s Manual that comes with the computer. You can always find something 
useful in it. 
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1. Resonant LC 
Circuit Parameters 


PROGRAM DESCRIPTION 


This program evaluates a resonant circuit having a lossless, single-layer airwound 
coil and a capacitor. The circuit is fully described by two equations—an approx- 
imation formula for the coil inductance and a resonance equation of an LC circuit. 


Formula 1-1 (2nF)? x LC = 1 p (№0)? x 10° 
I8D + 40A 


F is the frequency in Hertz, C is the capacitance in Farads, L is the coil inductance 
in Henrys, A is the coil length in inches, D is the coil diameter in inches, and N is 
the number of turns of the coil. 

The program helps in designing resonant circuits, winding coils, and repairing 
or modifying radio receivers and transceivers. When running the program first 
input all known parameters and then let the computer determine the missing ones 
by means of the above formulas. As shown in the examples, the program will solve 
fairly complicated problems. The parameters F, C, L, A, D, or N may beentered, 
computed, or reviewed in arbitrary order by pressing the appropriate keys. 


INSTRUCTIONS 


The program prompts with a question mark (?) for a command or with the name 
of the function when it is expecting a value. When prompted press one of the 
following keys: 
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S Clear all variables, start a new run 

F Enter frequency F in Hertz 

ZF Compute frequency from C and L 

C Enter capacitance C in Farads 

AC Compute capacitance from F and L 

L Enter coil inductance L in Henrys 

ZL Compute inductance from F and C 
ZX,orX Compute coil inductance from A, D, and N 
A Enter coil length A in inches 

ZA Compute coil length from L, D, and N 

D Enter coil diameter D in inches 

ZD Compute coil diamter from L, A, and N 

N Enter number of turns N for the coil 

ZN Compute number of turns from L, A, and D 


To review value of any parameter, key Z followed by F, C, L, A, D, or N. 


EXAMPLES 
Problem Definition 


1. Design a resonant circuit using a variable capacitor with capacitance in the 
range of 35-pF to 350-pF and a coil to be wound on a form 0.6-in. long and 0.3-in. 
in diameter. The lowest frequency should be approximately 6-MHz. Find the 
number of turns on the coil and the highest frequency at which the circuit will 
resonate with the capacitor at its minimum value. 

2. Asingle-layer coil 1.5-in. long, 0.6-in. in diameter, with 14 turns as shown in 
Figure 1-1 resonates with an unknown capacitor at 10-MHz. What parallel 
capacitor value is required to lower the resonant frequency to 9-MHz? 


14 
TURNS 


T | 1.5" Td 
Y 


Figure 1-1 
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Sample Run 

?F ФЕ 

FREQ.? FREQ.? 
FREQ.=6000000 FREQ.=10000000 
?C тА 

САР, ? LENGTH? 

CAP .=3.5E-10 LENGTH=1.5 

?7 тр 

СОМРОТЕ ШАМ, ? 

?L TIAM.=0.6 
INDUCT .=2.0103409E-6 ?N 

7A TURNS? 

LENGTH? TURNS=14 
LENGTH=0.6 ?X 

?П INDUCT .=9.9661017E-7 
LIAM. ? ?7 

DIAM.-0.3 COMPUTE 

TZ тС 

COMPUTE САР.-2.5416453Е-10 
TN ? 
TURNS=25.626381 

TN 5, 

TURNS? LET T=C 
TURNS=25 

7X ?F 

INDUCT .=1.9132653E—-6 ҒКЕО.? 

?7 ЕКЕй.=9000000 
COMFUTE TZ 

?7F COMFUTE 
FREQ.=6150331.3 тр 

"C САР.-3.1378337Е-10 
САР.? ? 

САР + =3 + 5Е-1 1 

ЕЕ РКІМТ С-Т 
ее 5.9518838Е-11 


FREQ.=19449055 


Discussion of Results 


The first problem yields 25.6 turns. In a practical application this number is 
rounded down to 25 turns. The recomputed resonant frequency is then 6.15-MHz. 
The highest frequency at which the circuit will resonate with a 35-pF capacitor is 
then found to be 19.44-MHz. 
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In the second problem we find a resonant capacitor of 254-pF at 10-MHz. We 
then stop program execution with BREAK and store this value temporarily as an 
arbitrary variable T. Program execution is then resumed with the CONT key and 
the resonating capacitance at 9-MHz is found (313-pF). The program is stopped 
again and the difference is found to be 59-pF, the final answer. This example illus- 
trates alternate use of the Timex/Sinclair 1000 in the programming and imme- 
diate (calculator) mode. 


PROGRAMMING REMARKS 


The program illustrates the use of the INKEY$ sequence starting in line 40. The 
CODE function in line 95 returns a number (see page 137-139 of the User’s 
Manual) for each key pressed during the INKEY$ sequence. The program flow is 
then directed to a line number corresponding to the key pressed (line 105). If the 
key is not recognized an error trapping routine returns a new prompt in line 30. A 
flag (Z) is set to 0 or 1 depending on whether the program expects an input (fre- 
quency, inductance, etc.) or whether it should find a missing parameter from the 
given ones. 


PROGRAM LISTING 


20 CLS 

22 LET Z=0 

25 SLOW 

27 PAUSE 20 

JO PRINT "т"; 

40 IF INKEY$-"" THEN GOTO 40 

45 PAUSE 10 

90 LET A$-INNEY$ 

60 PAUSE 5 

65 PRINT А% 

70 IF AS=INKEY$ THEN GOTO 70 

95 LET К-СОПЕ A$ 

100 IF K-56 THEN GOTO 20 

105 IF K=63 OR К-40 OR К-49 OR K=43 THEN 

GOTO KX10 
106 IF K=41 OR К-51 ОК K=38 ОК К-61 THEN 
GOTO KX10 

110 GOTO 30 

380 IF Z=1 THEN GOTO 384 

581 PRINT "LENGTH?" 

382 INPUT A 

383 GOTO 386 

384 LET 7-0 

385 LET A-CCOIXNX1E-30XX2-18XIxXL )/40/L 
384 PRINT "LENGTH-" sf 


387 
400 
401 
402 
403 
404 
405 
406 
407 
410 
411 
412 
413 
414 
415 


416 
417 
450 
431 
432 
453 
434 
435 
436 
437 
490 
491 
492 
493 
494 
495 
496 
497 
210 
Jil 
912 
213 
214 
515 
216 
517 
518 
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GOTO 30 

IF 7=1 THEN COTO 404 
PRINT "CAF, ?" 

INFUT C 

GOTO 406 

LET Z=0 

LET Cz1/C2XPIXF )kk2/L 
PRINT "САР. =" С 

GOTO 30 

IF 7-1 THEN GOTO 414 
PRINT "“DIAM.?" 

INPUT I 

GOTO 416 

LET Z=0 

LET I-(CT?XL4SQR ( 81 XLKXL+40KLKAKNEN/IE-6) 
)/N/NX1E6 

PRINT "HIAM,=" yl 
GOTO 30 

IF Z=1 THEN GOTO 434 
PRINT "FREQR, ?" 

INFUT F 

GOTO 435 

LET Z=0 

(ЕТ F=1/2/FI/SQR (LXC) 
PRINT "FREQ.=" $F 
GOTO 30 

IF Z=1 THEN GOTO 494 
PRINT “ІМІШСТ.?" 
INFUT L 

GOTO 496 

LET Z=0 

LET L-1/C2XPFIXFOXX2/C 
PRINT "INIUCT.z-"sL 
GOTO 30 

IF Z-1 THEN GOTO 514 
PRINT "TURNS?" 

INPUT N 

GOTO 517 

LET 2-0 

GOSUR 1000 

LET N-SQR (XXL )/I! 
PRINT "TURNS-z" 5N 
GOTO 30 


5 


6 


610 
611 
612 
613 
630 
651 
632 
1000 
1010 
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LET 7-0 

GOSUR 1000 

LET L-CNXIDDXX2/X 

GOTO 496 

LET 2-1 

PRINT "COMPUTE" 

GOTO 30 

LET Х-( 18XIt-40XAOX1E6 
RETURN 


Originally published as “Pocket Computer Solves for Resonance Using BASIC.” 
Reprinted from Electronics, June 16, 1981, copyright € McGraw-Hill, Inc. 1981. 
АП rights reserved. 


2. Resistive Attenuation 
and Matching Pads 


PROGRAM DESCRIPTION 


This program finds resistive component values for attenuation and matching pads 
commonly used in transmission networks. It computes series and shunt resistors 
for six types of resistive pads; unbalanced T and Pi, balanced O and H, and the 
minimum loss matching unbalanced L or balanced C pads. In addition to finding 
the resistor values, the program also determines the minimum loss of an imped- 
ance matching pad. If DC insulation is not required an impedance matching pad is 
often used as an inexpensive and simple substitute for an impedance matching 
transformer. Designation of the components is shown in Figures 2-1, 2-2, and 2-3. 
If the program finds negative resistor values for a requested input impedance, out- 
put impedance and loss, then that pad cannot be realized without gain elements 
and the requested loss should be increased. 


Figure 2-1 Unbalanced T/Balanced Н RI/2 R2/2 








RI R2 
R3 
Z | — ~<—72 
R3 
RI/2 R2/2 
N == 10-055 (dBy/10 R2 = 72 N + l 12% R3 
N-—1 
9 VNZ x Z2 
DENEN. DE RI = a - вз 
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Figure 2-2 Unbalanced тт, Balanced О 


R3/2 

















М = 1 21 ж 72 1 I IN + 1 l 
R3 = "cM NR Heat жи 
- N R2 Z2lN- 1] R3 
р Ере oA 
КІ ЖАДЫ — 1 R3 
Figure 2-3 Minimum Loss 
RI RI/2 
"АБ -- 21-->- R3-—Z22 
RI/2 


Е 72 Е 71 + 71 E 
КІ = ZI Vi-7: TEE V 2-1) 


ME ЧН. таа Loss (db) = 10 log М 


Ll 22 
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INSTRUCTIONS 


After pressing the RUN key the computer prompts for the type of attentuation/ 
matching pad. Provide the proper code, 1 for an unbalanced T or a balanced H, 2 
for unbalanced Pi or balanced O, or 3 for a minimum loss balanced C or unbalanced 
L matching pad between different impedances. The next prompts are for Z1 (in- 
put impedance), Z2 (output impedance), and loss, except for a minimum loss 
matching pad, where the loss is computed by the program. Units for all inputs and 
outputs are ohms and decibels. 


EXAMPLES 
Problem Definition 


1. Given input impedance Z1 = 600-ohm, output impedance Z2 = 900-ohm, 
find component values for a balanced Pi pad with a loss of 15-db. This type of pad 
could be encountered in a telephone office. 

2. Your TV is overloaded with too strong a signal. Design a balanced O pad (4 
resistors) for 71 = Z2 = 300-ohm and 10-db of loss. 

3. Design a minimum loss pad between 75 and 300-ohm impedances. You could 
use it as a match between a coaxial cable and the terminals of your TV set. 


Sample Run 


ATTENUATION PALS 

UNREAL. Ts BAL. PI - ENTER 1 
MIN. LOSS - ENTER 3571 

ENTER 717600 

ENTER 727900 

ENTER LOSS IN П6?15 
R1=369.29931 OHM 
К2-688.89257 OHM 
К5=269,.88721 OHM 


ATTENUATION PALS 

UNBAL. Т» BAL. PI ~ ENTER i 
UNBAL. РІ» RAL. Ü — ENTER 2 
MIN. LOSS - ENTER 372 

ENTER Z1i?300 

ENTER Z273500 

ENTER LOSS IN 06710 
R1=577.48518 OHM 
К2-577.48518 OHM 
R3=426.90748 OHM 
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ATTENUATION FADS 

UNBAL. Т, BAL. PI - ENTER 
UNBAL. РІ» BAL. 0 — ENTER 
MIN. LOSS - ENTER 373 

ZL > £2 

ENTER 71%300 

ENTER Z2?75 
LOSS=0.57194735 ПЕ 
R1-259.80762 OHM 

R3=86, 60254 OHM 


м н 


ATTENUATION FADS 

UNBAL. Т» BAL. PI - ENTER 
UNBAL. РІ, BAL. O — ENTER 
MIN. LOSS - ENTER 3? 


М t 


Discussion of Results 


In general, standard resistor values can be found within a few percent of the com- 
puted values. Notice that in the third example the minimum loss is only 0.57-db. 
This would compare with between 1 and 2-db for a matching transformer. 


PROGRAMMING REMARKS 


Error checks for the input values are performed in line 75. The proper formulas 
are then selected depending on the type of pad chosen. To save code the same in- 
put sequence (lines 90-125) for Z1 and Z2 is used for all three types of pads. 
Similarly the same output (print) sequence (lines 350-370) is used for all three 
branches. 


PROGRAM LISTING 


10 CLS 

30 PRINT "ATTENUATION PADS" 

40 FRINT "UNRBAL. T» BAL, PI — ENTER 1" 

50 PRINT "UNRAL. PI» BAL. 0 - ENTER 2" 

60 PRINT "MIN. LOSS - ENTER 5%"; 

70 INPUT А 

72 PRINT А 

75 IF A<>INT А OR А>3 OR A<1 THEN COTO 280 
80 IF A=3 THEN FRINT "Z1 > Z2" 

90 PRINT "ENTER 717"; 


100 
105 
110 
120 
125 
130 
140 
150 
155 
160 
170 
180 
190 
200 
210 
220 
230 
240 
290 
260 
270 
280 
290 
300 
310 
320 
540 
350 
360 
570 
375 
380 


INPUT 
PRINT 
PRINT 
INFUT 
PRINT 


IF A=3 THEN GOTO 270 
"ENTER LOSS IN DBT" $ 


PRINT 
INPUT 
PRINT 


LET N-10x*t( 0/10) 
IF А>1 THEN GOTO 
LET G=2*SQR ( МЖБЖС)/( М-1) 
LET Е-ВЖСМ412/(М-19-6 
LET F=CX( N+1 9/(N-10-G 


GOTO 
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B 
B 


"ENTER 72%"; 


C 
C 


I! 
n 


350 


220 


IF A=3 THEN COTO 270 


LET б-(М-12/2ж50К CEXC/N) 
LET E-1/CCNT1)D/CN-1)2/B-1/6) 
LET Fz1/CKCNT1)2/CN-1)0/C-1/6G) 


GOTO 


IF 1-С/В>0 THEN GOTO 300 
“ЕККОК» REENTER" 


PRINT 
GOTO 


550 


40 


LET EZEXSQR (1-С/В) 
LET G=C/SQR (1-С/В) 


LET L-LN (SRR (Б/СОФ5БОК (В/С-102020ЛМ 10 
"L0SS-"sLs" BB" 


PRINT 
PRINT 


IF А<>3 THEN PRINT "R2-"sFs" OHM" 


PRINT 
PRINT 
GOTO 


"RI="$E5" 


"R3z" $65" 


20 


OHM" 


OHM" 
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3. Characteristic Impedence 
of Transmission Lines 


PROGRAM DESCRIPTION 


This program finds the characteristic impedance for the most common types of 
transmission lines — the single conductor coaxial, the balanced shielded line, the 
open two-wire line in air, the parallel strip line used at microwave frequencies, and 
the transmission line consisting of two shielded parallel wires with sheath return. 
The five types of transmission lines are shown in the following five figures with the 
appropriate equations. Units of length are arbitrary (inches, feet, centimeters, 
etc.), as only ratios of linear measures are used in the equations. 


1. Single coaxial line 


Figure 3-1 
138 D 
ZO = — log |5 
| ve > (2 
R О e—relative diel. constant 


Z | 


2. Balanced shielded line 


_ 70 L — й 
S. | 20 = 28 ор (2v E 


Figure 3-2 
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3. Open two-wire line in air 


Figure 3-3 


gs "н + VE 
Los 


4. Parallel strip line 


Figure 3-4 


b 1 Т ae < 
IN 


5. Two wires in parallel with sheath return 


99) 
з 
- 
zm 


a 


Figure 3-5 Same Figure as Figure 3-2 


че Г Ë А ‘ 
70 = e ӘБ a (1 — 0) 


с = H/D У = Н/К 
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INSTRUCTIONS 


Upon pressing the RUN key the program prompts for the type of transmission 
line; answer with 1, 2, 3, 4, or 5 referring to Figures 3-1 through 3-5. The next 
prompt is for the relative dielectric constant of the medium in which the conductor 
is embedded. This constant is 1.0 for air and is normally supplied by the dielectric 
manufacturer. Next, the program prompts for the diameter of the inner conduc- 
tor, spacing between conductors for a two-wire line, and shield diameter when ap- 
propriate. Refer to the above figures when entering values. When all inputs are 
given the program computes the characteristic impedance of the line in ohms. 


EXAMPLES 


Problem Definition 


1. Find characteristic impedance 70 of a single coaxial line with an inner con- 
ductor diameter of 1.5-mm and shield diameter of 6.5-mm. Assume a relative 
dielectric constant of 1.4 in this and the following examples. This would be a 
typical coaxial cable used for TV downlead from the antenna. 

2. Find ZO for a balanced shielded cable with conductor spacing of 10-mm, 
outer diameter of 14-mm and conductor diameter of 1.2-mm. 

3. Find ZO for a two wire line with a spacing of 9-mm and a conductor diameter 
of 1.5-mm. This would be a typical twin lead cable used for TV. 

4. Find ZO of a parallel strip line, of the type used in microwave transmission, 
with conductor width of 11-mm and spacing of 1-mm. 

5. Finally, find ZO of two parallel lines with sheath return. The spacing is 
10-mm, the outer diameter is 14-mm and the inner diameter is 1-mm. 


Sample Run 


CHAR. ІМРЕПАМСЕ 


ENTER ТҮРЕ {1-5 )?1 

ENTER REL. DIAL. CONSTANT?1.4 
SINGLE COAX. LINE 

OUTER DIAM П?6.5 

COND. DIAM. RPI.S 
20-74.275581 OHM 
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ENTER TYPE (1-5 )?2 

ENTER REL. DIAL. CONSTANTT1.4 
BAL. SHIELDED LINE 

SPACING Н?10 

OUTER DIAM 1714 

COND. ШАМ. RT7T1.2 
20-170.94141 OHM 


ENTER ТҮРЕ (1-53 

ENTER REL. DIAL. CONSTANTT1.4 
2 WIRE LINE 

SPACING H?9 

CONI, DIAM. К?1.5 
20-251.02165 OHM 


ENTER TYPE (1-5274 

ENTER REL. DIAL. CONSTANT?1.4 
PARALLEL STRIP LINE 

CONDUCTOR WIDTH W?11 

SPACING Н?1 

20-28.965741 OHM 


ENTER ТҮРЕ {1=9 779 

ENTER REL. DIAL. CONSTANT?1.4 
2 FAR. LINES» SHEATH RET. 
SPACING Н?10 

OUTER DIAM IT714 

CONI. DIAM. RPI 

20-50.1676 OHM 


Discussion of Results 


As shown in the listing of the run the characteristic impedances for the five ex- 
amples are 74.3, 170.9, 251.0, 29.9, and 50.2-ohm. 


PROGRAMMING REMARKS 
Note error checks in lines 37, 120, 220, 320, 450, and 520. The same subroutine 


starting in line 600 is used for input to various branches of the program to save 
steps. 
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PROGRAM LISTING 


15 CLS 
20 PRINT "CHAR. IMPEDANCE" 
25 PRINT 
27 PRINT "ENTER TYPE (1-5)7"; 
30 INPUT A 
35 PRINT А 
37 IF А<2ІМТ А ОК А<1 OR А>5 THEN GOTO 60 
40 PRINT "ENTER REL. DIAL. СОМ5ТАМТТ"; 
45 INPUT E 
30 PRINT E 
Jo IF E20 THEN GOTO 100жА 
60 LET L=25 
65 GOTO 720 
100 FRINT "SINGLE COAX. LINE" 
110 GO0SUR 650 
120 IF К>0 AND DSR THEN GOTO 150 
130 LET L=110 
140 GOTO 720 
150 LET Z=138/SQR ЕЖ М (CI/R)O/LN 10 
160 GOTO 380 
200 PRINT "BAL. SHIELDED LINE" 
210 GOSUR 600 
220 IF R20 AND H-R AND D2HtR THEN GOTO 250 
230 LET L=210 
240 GOTO 720 
250 LET V=H/R 
260 LET S=H/D 
270 LET 7=276/5ЧК EXLN (2XUX( 1-5Ж52/ 
(1+SkS))/LN 10 
280 GOTO 580 
500 PRINT "2 WIRE LINE" 
310 GOSUR 600 
320 IF R20 AND НУК THEN GOTO 350 
330 LET L=310 
340 GOTO 720 
350 LET S=H/R 
360 LET Zz276XLN (StSQK ( 5Ж5-102/15 10/5QK E 
370 GOTO 580 
400 PRINT "PARALLEL STRIF LINE" 
410 PRINT "CONDUCTOR WIDTH Ыт"; 
420 INPUT W 


430 
440 
450 
460 
470 
480 
490 
300 
210 
320 
530 
340 
550 
560 
370 
580 
290 
600 
610 
620 
630 
640 
650 
640 
670 
680 
690 
700 
710 
720 
730 
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FRINT W 

GOSUBR 600 

IF Н>0 АМП W-O THEN GOTO 480 
LET L=410 

GOTO 720 

LET Z=377XH/W/SQK E 

GOTO 580 

PRINT "2 РАК. LINES» SHEATH RET." 
GOSUR 600 

IF R>O AND НЕК ANI П>Н+К THEN COTO 550 
LET L-510 

GOTO 720 

LET V=H/R 

LET 5-Н/П 

LET Z-69/5GR EXLN (U/2/5/S5X(1-SX*40)0/LN 10 
PRINT "ZO="5Z35" OHM" 

GOTO 25 

FRINT "SPACING H?"; 

INPUT H 

FRINT H 

IF 4-3 THEN GOTO 680 

IF 4-4 THEN RETURN 

PRINT “OUTER DIAM LD?" $ 

INFUT [I 

PRINT D 

PRINT "COND. DIAM. R?"; 
INPUT К 

PRINT К 

RETURN 

PRINT “ERROR: КЕЕМТЕК" 

GOTO L 





ОАТА 
TRANSMISSION 


4. Bit Error Rates for Various 
Modulation Schemes 





PROGRAM DESCRIPTION 


In a well-designed data transmission link, random noise with Gaussian distribu- 
tion is the main source of data bit errors. The probability of receiving an incorrect 
data bit depends on both the Signal-to-Noise (S/N) ratio and the type of data 
modulation used. In general, circuit complexity is traded for error probability. A 
more complex circuit with tighter attenuation and delay distortion requirements 
will result in a lower probability of error for the same S/N ratio. The following 
program computes the Bit Error Probability (Pe) as a function of S/N for nine 
popular modulation schemes. The formulas used are as follows: 


1. On-Off Keying (OOK), Coherent: 


Formula 4-1 Pe 


2. Frequency Shift Keying (FSK), Coherent 
Amplitude Shift Keying (ASK), Coherent 
Pulse Code Modulation (PCM), Unipolar 


E 
Formula 4-2 Pe = 5 Erfc 


- 
Zu 
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3. PCM, Polar 4. FSK, Noncoherent 
Phase Shift Keying (PSK) 


Formula 4-3 Formula 4-4 
5d 5 214%) 
Ре "t у, Erfc N Pe = 2 е 


5. Differential Phase Shift Keying (DPSK) 6. ASK, Noncoherent 


Formula 4-5 Formula 4-6 


The function Erfc (complementary error function) related to the Gaussian noise 
distribution is defined as follows. The function is computed recursively, as shown 
in Program 12. 


Formula 4-7 Erf (x) = = | e dt Erfc (x) = 1 — Erf (x) 
0 


INSTRUCTIONS 
When the RUN key is pressed the prograt prompts for the signal-to-noise ratio in 
db and the type of modulation. Select a S/N ratio and a number between 1 and 6 


to indicate the type of modulation as described in the Program Description. The 
program then computes the bit error probability. 


EXAMPLES 


Problem Definition 


Find the probability of bit error Pe for 12-db S/N for all nine modulation 
schemes. 
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Sample Run 


BIT ERROR PROBABILITY 

ENTER S/N RATIO ІМ ПЕ?12 
SELECT MODULATION CODE 1-671 
РКОВ( 12 ПЕ )=.0024385373 


ENTER S/N RATIO IN ПЕ?12 
SELECT MODULATION CODE 1-672 
РКОВ( 12 ПВ )=.000034 302624 


ENTER S/N RATIO IN ПЕ?12 
SELECT MODULATION CODE 1-473 
PROBC12 ВЕ2-9.0060101Е-9 


ENTER S/N RATIO IN IRT12 
SELECT MODULATION CODE 1-674 
РКОВ( 12 IB2-,00018089715 


ENTER S/N RATIO IN ПЕ?12 
SELECT MODULATION COME 1-675 
РКОВЕ( 12 IB0-6.,544347E-8 


ENTER S/N RATIO IN ПЕ?12 


SELECT MODULATION CODE 1-676 
PROE 12 DB)=,.0001970186 


Discussion of Results 


The probability of error for the nine modulation schemes varies widely. The 
lowest probability of error is given by type 3 (PCM, Polar/PSK) modulation 
scheme, the highest by type 1 (OOK, Coherent) modulation scheme. 


PROGRAMMING REMARKS 


The seemingly different formulas are solved by the same section of code through 
recomputing the argument of the complementary error function. Further discus- 
sion of recursive evaluation of the error function is given in Program 12. 


PROGRAM LISTING 


15 CLS 


20 PRINT "ВІТ ERROR PROBABILITY" 
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PRINT "ENTER S/N RATIO IN IIR?" ; 
INPUT B 

PRINT В 

LET C=B 

PRINT "SELECT MODULATION COIE 1-67"; 
INFUT D 

PRINT D 

IF П>0 AND 047 AND D-INT D THEN GOTO 100ЖП 
PRINT "ENTRY ERROR" 

GOTO 45 

GOSUR 800 

GOSUR 800 

LET Х-10жж( C/20) 

LET Yz1/SüQR РІЖЕХР (-XXX) 

IF Х<2 THEN GOTO 350 

LET А=14/Х 

FUR 1527 TU 1 STEP -1 

LET A-I/2/CXtà) 

NEXT I 

LET Z=Y/( X+ñ ) 

LET М=1-7 

GOTO 700 

LET I=1 

LET Т=2ЖХЖҮ 

LET N=T 

LET Т-2ЖХЖХЖТ/(2Ж141) 

LET М-МҒТ 

IF M=N THEN GOTO 393 

LET N=M 

LET I=I+1 

GOTO 365 

LET Z-1-N 

GOTO 700 

GOSUR 800 

LET Z=EXP -C10XXC(C/10)) 

GOTO 700 

LET F-10X*(C/10) 

LET Z-EXF (-Е/2)Ж( 1+1/5ЧК C2XPIXF 5) 


LET Z=Z/2 

PRINT "РКОВС"%В;" UDB)=" +Z 
РКІМТ 

GOTO 25 

LET C-C-10XLN 2/LN 10 


RETURN 
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NUMBER THEORY 


5. Any Base fo Any Base 
Conversion 


PROGRAM DESCRIPTION 


Converting numbers between different bases is frequently done by computer pro- 
grammers, mathematicians, high school, and college students. In particular, com- 
puter program development requires conversion between binary (2), octal (8), 
decimal (10), and hexadecimal (16) bases. Number theory problems may require 
conversion between these and other bases. 

Conversion between two arbitrary bases a and b is generally a two-step pro- 
cedure, converting first from base a to base 10 and then from base 10 to base b. 
The unique method presented here for positive integers uses the same algorithm 
for both conversions, thus saving on programming steps. 

To represent a number in a base larger than 10, more than one place is necessary 
to represent one digit in that base. For example 17FE (base 16) would be entered or 
displayed as [01], [07], [15], [14] or simply 1071514 (1=01, 2=02,..., A=10, 
B=11,C=12, D=13, E=14, F= 15). In general, one place would be reserved per 
digit for bases 2-10, two places for bases 11-100, etc. The program automatically 
interprets the entry or display depending on the base. The program has no in- 
herent limitations on the magnitude of the base or the number to be converted. 
The only limiting factor is the computer’s accuracy. For example, the accuracy of 
the Timex/Sinclair BASIC will limit the input and output to about 10 decimal 
places. A larger number would result in error. 


INSTRUCTIONS 


After pressing the RUN key follow the prompts by entering the number to be con- 
verted (argument), and the old and the new base. The program will display the 
argument converted to the new base; if neither the old nor the new base is equal to 
10, then the argument in base 10 will also be displayed. 


EXAMPLES 


1. Convert 17C (base 16) to base 2. 
2. Ronvert the decimal equivalent of 17C to base 2 (both results should give the 
same answer). 
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3. Convert decimal argument 12345 to base 122. 
4. Convert it back from base 122 to base 10. 


Sample Run 


GENERAL BASE CONVERSION ENTER ARGUMENTT12345 


ENTER ARGUMENT?10712 ENTER OLI BASE?TIO 
ENTER ШШ BASE?16 ENTER NEW RBA5ET122 
ENTER NEW RASET? 12345 TO BASE 10 
10712 TO BASE 16 =101023 TO BASE 122 
=101111100 TO BASE 2 
=380 TO BASE 10 ENTER ARGUMENT ?P1I01023 
ENTER QLI! BASE?TI22 
ENTER АКСОМЕМТ? 380 ENTER NEW EFASE?16 
ENTER OLIE BASE?10 101023 TO BASE 122 
ENTER NEW BASE?2 =3000307 TO BASE 16 
380 TO BASE 10 -12345 TO BASE 10 


=101111100 TO BASE 2 


Discussion of Results 


Some interpretation of results is required for bases larger than 10. The result of the 
third example should be interpreted as [101], [023] (base 122). 


PROGRAMMING REMARKS 


Extensive error checking of input data is performed in lines 70, 110, and 150. 
Round-off error is minimized by using the INT function and adding 0.5 to the ap- 
propriate functions in lines 440 and 480. Line 480 finds the number of decimal 
digits required to represent each ‘‘digit’’ to an arbitrary base. 


PROGRAM LISTING 


20 CLS 

30 PRINT "GENERAL BASE CONVERSION” 
40 PRINT "ENTER ARGUMENT?" ; 

90 INFUT С 

60 PRINT С 

70 IF C<>INT (ABS С) THEN GOTO 40 
80 PRINT "ENTER DOLI! BASE?" 3 

70 INFUT D 
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100 PRINT D 

110 IF D<>INT (ABS It) THEN GOTO 80 

120 PRINT "ENTER NEW BASE?"; 

130 INPUT E 

140 PRINT E 

150 IF E<>INT (AES E) THEN GOTO 120 

160 IF П<>10 THEN GOTO 190 

170 LET N-C 

180 GOTO 250 

190 LET U-I 

200 GOSUR 480 

210 LET S-C 

220 LET Q=D 

230 LET R=V 

240 GOSUR 390 

250 IF E<>10 THEN GOTO 280 

260 LET S=N 

270 GOTO 340 

280 LET U=E 

290 GOSUR 480 

300 LET S-N 

310 LET Q=V 

320 LET R-E 

330 GOSUR 390 

340 PRINT C$" TO BASE "ПІ 

350 РКІМТ "-";N$" TO BASE "; 

360 IF П<>210 ANI Е<>10 THEN PRINT "="$53" 
TO BASE 10" 

370 FRINT " " 

380 GOTO 4 

390 LET M=0 

400 LET М-0 

410 LET Ғ-5 

420 LET T-P 

430 LET F-INT (P/R) 

440 LET N=N+INT ((T-PXROXQXXM4,5) 

450 IF Р=0 THEN RETURN 

460 LET М=М+1 

470 СОТО 420 

480 LET U-INT (10XX( 1+INT (LN CU-1)/LN 10))+,5) 

490 RETURN 


6. Hex-to-Decimal and 
Decimal-to-Hex Conversion 


PROGRAM DESCRIPTION 


This program differs from the previous all-base conversion program in that it con- 
verts only between bases 10 and 16. The two conversions (10 to 16 and 16 to 10) are 
very important for computer program development as most programming tools 
(assemblers, memory maps, disassemblers, etc.) refer to memory locations and 
memory contents in either base 10 or 16 (hex). This program, unlike the previous 
one, uses the standard alphanumeric presentation for hexadecimal digits: 0 
through 9 and A, B, C, D, E, and F for digits 10 through 15. The same comments 
about accuracy and round-off as in the previous program apply. 


INSTRUCTIONS 


Press the RUN key and follow prompts to perform Hex-to-Decimal (1), or 
Decimal-to-Hex (2) conversion. After the next prompt enter the argument either 
in Hex (using A through F if necessary) or as a decimal number. 


EXAMPLES 


1. Convert F697 (base 16) to base 10. 
2. Convert the result back to base 16. 
3. Convert FFFF (base 16) to base 10. 


Sample Run 


НЕХ-Т0-ПЕСІНЕА! + ENTER 1 
DECTMAL-TO-HEXs ENTER 271 
ENTER ARGUMENT ІМ HEXTF697 
=63127 IN BASE 10 


HEX-TO-DECTHAL + ENTER 1 
QNECIMAL-TO-HEXs ENTER 272 

ENTER ARGUMENT TO BASE 10765127 
-Ғ697 IN BASE 16 
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HEX-TO-DECTMAL + ENTER 1 
DECIMAL-TO-HEXs ENTER 271 
ENTER ARGUMENT IN HEXTFFFF 
=45535 IN BASE 10 


Discussion of Results 


Prompting and results are self evident. 


PROGRAMMING REMARKS 


Error checking of input data is performed in lines 70 and 270. Base 16 digits are 
entered as a character string (A$) and each character in the string is interpreted in 
line 260. If a character is found which does not correspond to a number between 0 
and 15 it is rejected as an error. Translation from decimal numbers into a hex 
character string occurs in line 510. This program demonstrates the power of 
BASIC character string manipulation functions such as CHR$, CODE, and LEN. 


PROGRAM LISTING 


10 
20 
350 
35 
40 
30 
60 
70 
80 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
3510 


DIM В%(10) 

CLS 

FRINT 

FRINT "HEX-TO-DECIMALs ENTER 1" 


PRINT "IIECIMAL-TO-HEX» ENTER 27%; 


INPUT Е 

PRINT В 

IF B<>1 AND Е<>2 THEN GOTO 30 
GOTO 200ЖЕ 

PRINT "ENTER ARGUMENT IN HEX?" 
INPUT АФ 

PRINT АФ 

LET L-LEN f$ 

LET N=0 

FOR І-0 TO L-1 

LET С=СОПЕ A$L-I2-28 

IF С>=0 ANI С<=16 THEN GOTO 310 
FRINT "БАП ENTRY" 

PAUSE 200 

RUN 

LET М-МҒСЖ16ЖЖТІ 


320 
530 
340 
400 
410 
420 
430 
440 
450 
460 
470 
480 
485 
490 
500 
510 
380 
290 
600 
610 
620 
630 
640 
650 
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NEXT I 

PRINT "-"jN$" IN BASE 10" 
GOTO 640 

PRINT "ENTER ARGUMENT TO BASE 107"; 
INPUT D 

PRINT В 

IF П-ІМТ (ABS D) THEN GOTO 470 
PRINT "BAD ENTRY" 

PAUSE 200 

RUN 

LET К=0 

LET К=К+1 

LET F=D 

LET D-INT (1/16) 

LET E-F-16XD 

LET B$(K)-CHR$ (E428) 

IF Е>0 THEN GOTO 480 
PRINT "-"j 

FOR I-1 TO K-1 

PRINT EB$(K-I); 

NEXT I 

PRINT " IN BASE 16" 

PRINT * " 

GOTO 30 
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7. Prime Factor 
Decomposition 


PROGRAM DESCRIPTION 


One of the first problems facing a student in number theory is the prime factor de- 
composition of an integer. Prime number decomposition is not merely a theoreti- 
cal exercise—it is applicable in many areas, including cryptography. Some of the 
modern cryptographic codes are based on the assumption that it is extremely time 
consuming, even for large computers, to find prime factors of large numbers. In 
the so-called **public Кеу” technique the encryption code (product of two primes) 
can be publicized, while the prime factors being the decryption code, are known 
only to authorized recipients. If a code is based on a key that is a product of two 
large primes, which it would take years for a large computer to find, then the code 
is for all practical purposes unbreakable. 

The two principal methods for finding prime factors of an integer are the 
Eratosthenes sieve, named after a famous Greek mathematician, and a simple odd 
divisor check. The first method is faster but also more appropriate for large com- 
puters, as it stores in memory all successive primes smaller than the square root of 
the number to be decomposed. The second method is much simpler, as it only 
checks for decomposition into 2, 3, 5, and the successive odd integers. Though the 
second approach is somewhat wasteful on execution time, since it checks against 
all odd but not necessarily prime integers, the second method imposes low 
memory requirements on the computer, because no intermediate results have to 
be stored. The second method is therefore the only one suitable for calculating 
primes on the Timex/Sinclair 1000 computer, and is the basis for the following 
program. 


INSTRUCTIONS 
Press the RUN key to start the program, then follow the prompt by entering the 


number to be decomposed into prime factors. The program then displays how 
many prime factors were found and the prime factors, one by one. 


EXAMPLES 


Problem Definition 


Decompose 12345, 66759, and 71 into prime factors. 
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Sample Run 

PRIME FACTOR FINDER PRIME FACTOR ҒІМПЕК 
ENTER NUHRER?12345 ENTER NUMBER?TSS759 
NO. OF FACTORS=3 NO. OF FACTORS=5 
FACTOR NO. 1=3 FACTOR NO. 1=3 
FACTOR NO. 2-5 FACTOR NO. 2-7 
FACTOR NO. 3-82 FACTOR NO. 3-11 


FACTOR NO. 4-17 
FACTOR М0. 5-17 


РКІМЕ ҒАСТОК ҒІМПЕК 
ENTER МИМБЕК?71 

NO. ОҒ FACTORS=1 
FACTOR М0. 1-71 


Discussion of Results 


Notice in the second example that 17 appears twice, indicating that it is a multiple 
factor. In the third example we test a prime number (71) which has a single non- 
trivial factor, itself. 


PROGRAMMING REMARKS 


The program performs error checking in line 90. If the number to be decomposed 
is larger than 1,000,000,000 then round-off errors will not allow the program to 
operate correctly. For user recall the first 50 prime factors are stored in a dimen- 
sioned array B(50) and could be recalled by a user written program. If there are 
more than 50 factors then the program displays only the first 50. 


PROGRAM LISTING 


20 CLS 

25 LET K=0 

27 LET I-0 

50 DIM В(50) 

40 PRINT "PRIME FACTOR FINDER" 
30 FRINT "ENTER NUMBER?" y 

60 INPUT А 

70 FRINT А 

80 LET F=A 


30 


90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
270 
500 
510 
512 
515 
520 
339 
540 
550 
360 
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IF A=INT A ANH А>1 ANI! А<1Е? THEN. GOTO 120 
PRINT "ЕККОК»КЕЕМТЕК" 

GOTO 50 

LET E=A/2 

LET С=2 

IF Е<>ІМТ E ОК 2*E<A THEN GOTO 170 
GOSUR 320 

GOTO 120 

LET П=П+1 

LET C=1+2xk0 

IF CSR A THEN GOTO 250 

LET E=A/C 

IF Ex>INT E OR ЕЖГ<>А THEN СОТО 170 
LET I-I-1 

GOSUR 320 

GOTO 170 

IF А<=1 THEN GOTO 280 

LET К=К+1 

LET BC KO=A 

PRINT " NO. OF FACTORS="sk 

FOR I-1 TO К 

PRINT "FACTOR NO. "sIs"z"sRECI) 
NEXT I 

PRINT * " 

GOTO 25 

IF К>=50 THEN RETURN 

LET К=К+1 

LET R(K)=C 

LET А=Е 

RETURN 





COMPUTER 
PROGRAMMING 





8. Relative Offset 
Computations for Machine 
Code Programming 


PROGRAM DESCRIPTION 


The purpose of this program is to compute relative addresses of machine code in- 
structions. When writing short programs in machine code for any of the popular 
8-bit microprocessors such as Z80, 6502, or 6800 one either has to load the 
Assembler program, which may not be available, or ‘‘hand’’ assemble the 
machine code program. The main problem then is to find the arguments of the 
relative branching instructions. The program shown here will find in hexadecimal 
notation the argument of a relative branching instruction zz given the absolute ad- 
dress of that instruction, the ‘‘FROM?’’ address (X) and the hexadecimal address 
of the instruction to which the program should branch if the branching condition 
is true, the “ТО” address (Y). For example, given the following code in 780 
machine language: 








Address Instruction Label Mnemonic 
20FC DD 7E 00 LOOP LD A,(X) 
2109 B8 CPB 
210A 20 zz JR NZLOOP 


one has to find the value of zz, also known as offset, given the address of the 
“FROM” branching instruction 210A (hex) and the address of the “ТО” instruc- 
tion 20FC (hex). To compute zz one of the following equations must be solved: 
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zz = Y-X-2for Y > X (forward branch), or 
zz = FF (hex) + Y—X-—1 for Y < = X (reverse branch). 


The resulting offset in the above example would be FO (hex) and the last instruc- 
tion would read 20 FO. Because the offset zz can only assume values between 0 and 
FF (hex) for 8-bit microprocessors, the relative branching is limited to one-half 
page (1 page = 256 locations) in the forward or reverse direction. The program 
checks whether the branching instruction can be performed, as it analyzes the ad- 
dresses of the FROM/TO instructions. The program then indicates if the relative 
branching can or cannot be performed in a single step, and if it will result in a for- 
ward or a reverse jump. 


INSTRUCTIONS 


When the RUN key is pressed the program prompts for the address of the 
“FROM” and the “ТО” instructions. Both should be entered in hexadecimal 
notation. 


EXAMPLES 


Problem Definition 


Branch Instruction Destination Instruction 
Example # Absolute address (hex) Absolute address (hex) 
1 123A 127B 
2 25F3 2613 
3 010A OOFC 
4 AA25 BA24 


Sample Run 


REL. ADDRESSING (8-BIT) 


ENTER "FROM" ADDRESS IN HEX7?1234 
ENTER "TO" ADDRESS IN HEX7T127E 
FORWARD OFFSET-S3F 


ENTER "FROM" ADDRESS IN HEX"?25SF 5 
ENTER "TO" ADDRESS IN НЕХ?2613 
FORWARD OFFSET=1E 
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ENTER "FROM" ADDRESS IN HEX?10A 
ENTER "TO" ADDRESS IN HEXTFC 
REVERSE OR IN FLACE OFFSET=FO 


ENTER "FROM" ADDRESS IN НЕХ?АА25 
ENTER "TO" ABIIRESS IN НЕХ?ВА24 
4095(< ПЕС) OFFSET TOO LARGEsREENT 
ER 


ENTER "FROM" АШПКЕ55 IN HEXTAA25 
ENTER "TO" ADDRESS IN НЕХ?АА23 
REVERSE OR IN FLACE OFFSET=FC 


Discussion of Results 


The first three examples result in hexadecimal offsets of ЗЕ, 1E and FO. The 
fourth example indicates that the offset is too large. We realize then that the 
“ТО” address should be AA23 instead of BA24 and end up with an offset of FC. 


PROGRAMMING REMARKS 


Notice the use of double quotes (Shift Q) in lines 45 and 120. To perform arith- 
metic on hexadecimal numbers the numbers are first converted to base 10 and the 
results are converted back to base 16. Conversion from decimal to hex and vice 
versa is performed in the same way as in program £6 by means of string manipula- 
tion functions. 


PROGRAM LISTING 


20 CLS 

30 PRINT "REL. ADDRESSING (8-БІТ)" 

40 PRINT 

45 PRINT "ENTER ""FROM"" ADDRESS IN HEX?"; 
90 INPUT C$ 

60 PRINT C$ 

70 GOSUB 330 

80 IF F=0 THEN GOTO 110 

90 PRINT "ERROR? REENTER" 

100 GOTO 40 
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110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
245 
250 
260 
265 
270 
280 
290 
500 
320 
550 
340 
350 
560 
370 
380 
590 
400 
410 
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LET J-X 

PRINT "ENTER ""TO"" ADIKESS IN HEX?"; 
INPUT C$ 

PRINT C$ 

GOSUR 330 

IF F=0 THEN СОТО 190 

PRINT "ERROR? REENTER" 

GOTO 120 

LET I-X-J 

IF П<-1297 ANI П>--126 THEN GOTO 230 
PRINT Ds"( ПЕС) OFFSET TOO LARGEsREENTER" 
GOTO 40 

IF П<2 THEN GOTO 260 

LET ZzI-2 

LET H$-"FORWARI" 

GOTO 270 

LET 7=254+П 

LET H$-"REUERSE OR IN PLACE" 
(ЕТ G=INT (2/16) 

LET E$-CHR$ (6428) 

LET М%-СНК% (CZ-16X*X6t28) 
РКІМТ H$s" OFFSET="sE¢sW$ 
GOTO 40 

LET F=0 

LET Х=0 

LET L-LEN C$ 

FOR I=0 TO L-1 

LET К-СОПЕ C$CL-I2)-28 

IF К<0 OR К>15 THEN LET F=1i 
LET X-zXtRX16XXI 

NEXT I 

RETURN 


COMPUTER 
GENERATED 
PLOTTING 





9. Curve Smoothing and 
Interpolation / Extrapolation 
with Lagrange Polynomials 


PROGRAM DESCRIPTION 


This program is the mathematical equivalent of a ‘‘French Curve.” Its main ap- 
plication is interpolation and extrapolation of measured results. Given Х(1), Ү(1), 
i=1,...,n, the program finds a unique polynomial of degree (n — 1) fitting these n 
points, and evaluates the polynomial for arbitrary values of the independent 
variable X. The approach differs from a least squares approximation in that the 
resulting polynomial fits the points X(i), Y(i) exactly. 

The following formula, attributed to the French mathematician Lagrange, 
describes the polynomial: 


Formula 9-1 n 
т X — x) 
s=l,s #i 


y (х) = 2, y, X 


2 


Il 


т Кр. — X.) 


5-1,6 #1 


35 


36 SCIENCE AND ENGINEERING PROGRAMS 


For practical applications it is advisable to use care in selecting only adjacent 
points to compute the polynomial so that wildly oscillatory curves are not ob- 
tained as extrapolations and interpolations of relatively smooth functions. In gen- 
eral, a cubic polynomial based on four nearby points (n= 4) will be satisfactory. 

An important feature of the program is that a moving interpolation ‘‘window’”’ 
can be established so that only k measured values (К < = n) closest to the inter- 
polation point will be included in the computation. 


INSTRUCTIONS 


Pressing the RUN key initializes the program and issues prompts for the number 
of points (n < = 20) and the successive X(i) and Ү(1) values of points in increasing 
order of X; the number of points in the **window"' and the value of X at which the 
function should be interpolated or extrapolated. The next option is to recompute 
the function for another X, change the ‘‘window’’ width, or start a new run. 


EXAMPLE 


Problem Definition 


As shown in Figure 9-1 six points have been measured. Assume a ‘‘window’’ 
width of four and evaluate the function at X = 2.1 (interpolation) and at X = 7.0 
(extrapolation). Then change the ‘‘window’’ width to five points and recompute 
the same two points. 


Sample Run 


LAGRANGE INTERPOLATION 
NO. PFOINTS?6 

ENTER Х(127-0.5 

ENTER Y(1)57?0 

ENTER Х(2270.5 

ENTER Ү(2270.7 

ENTER X(3)72 

ENTER Ү‹З)?2 

ENTER Х( 4973 

ENTER Ү(4272.6 

ENTER Х( 5274.5 

ENTER Ү(52%2.8 

ENTER Х(62%6 

ENTER Ү(6273.3 

NO. POINTS ІМ WINTOQW?T4 
FINI! Y(X); Х-?2,3 

Y( 2,3)=2, 20796 


CURVE SMOOTHING AND INTERPOLATION/EXTRAPOLATION 


NEW X isNEW WINDOW 2sNEW RUN 371 
FINI Y(X); X=?7 
YC7 )=4 , 4333333 


NEW X i NEW WINDOW 2 ,NEU RUN 372 
NO. POINTS IN WINDOW?S 
FIND Y(X); Х-?2.3 


NEW X 1,МЕШ WINDOW 2sNEW RUN 571 
FIND Ү(Х)» Х-?/ 
Ү(79)3-5.1909091 


NEU X 1>МЕМ WINDOW 2»NEU RUN 371 
FIND ҮСХ), Х-?4.5 
Ү 4.5 2072.8 


Figure 9-1 
LEGEND: 


X - MEASURED 
© - COMPUTED 


y (4 POINTS IN WINDOW) 
6 O- COMPUTED 
(5 POINTS IN a 
5 WINDOW) 
© 





Discussion of Results 
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The example shows that ‘‘window’’ width strongly affects extrapolated values. 
The interpolated value changed only from 2.20. . to 2.22. ., whereas the extrapo- 
lated value changed from 4.43. . to 5.19... The user has to decide which ‘‘win- 
dow’’ width is more appropriate. The example also shows that the Lagrange poly- 
nomial goes exactly through each given point. When Y(4.5) is computed by the 
program, it generates 2.8, which is the exact value originally entered for Y(4.5). 
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PROGRAMMING REMARKS 


The program checks for the correct entries in lines 80, 170, and 280. Branching is 
performed in three simple statements in lines 560-580. The Lagrange formula 
evaluation takes place in a double loop in lines 430-500. 


PROGRAM LISTING 


ә CLS 
20 [TIM X20) 
30 DIM Y(20) 
40 PRINT "LAGRANGE INTERFOLATION" 
30 PRINT "NO. POINTS?"s 
60 INPUT А 
70 PRINT А 
80 IF A<=20 ANI! 422 ANI A=INT ABS A THEN 
GOTO 110 
90 PRINT "ERROR" 
100 GOTO 50 
110 FOR I-1 TO A 
120 PRINT "ENTER X("sIs")?"; 
150 INPUT X 
140 PRINT X 
150 LET X(IO-X 
160 IF Ісі THEN GOTO 200 
170 IF X(IO-XCI-1)5 THEN GOTO 200 
180 PRINT "XCITi0oz-XCI)" | 
190 GOTO 120 
200 PRINT "ENTER Y("sIs")?"; 
210 INFUT Y 
220 PRINT Y 
250 LET Y(I)=Y 
240 NEXT I 
250 PRINT "NO. POINTS IN WINDOWT" y 
260 INFUT R 
270 FRINT B 
280 IF B<=A ANI Е>і ANI B-INT ABS E THEN 
GOTO 310 
290 PRINT "ERROR" 
500 GOTO 250 
510 LET C-INT (R/2) 
520 LET I-1 
330 PRINT "FIND Y(X); X-?"; 
540 INFUT W 


350 
350 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
540 
570 
580 
590 


CURVE SMOOTHING AND INTERPOLATION/EXTRAPOLATION 


PRINT W 

IF W“<=X( I) OR I=A THEN GOTO 390 
LET I=I+1 

GOTO 360 

LET 2=0 

LET Ш=1-С 

IF П<1 THEN LET D-1 

IF I~CtB>A THEN LET П=А-В+1 
FOR Т=П TO П+Е-1 

LET Е=1 

FOR J=I TO D+E-1 

IF J=I THEN GOTO 480 

LET E-EXCU-XCJ))9/CXCIO-XCJI 2) 
NEXT J 

LET Z-ZFTEXYCI) 

NEXT I 

PRINT "Y( "yWy" D="5Z 

PRINT 


РКІМТ "NEU X 1,МЕШ WINDOW 2»NEU RUN 37"; 


INFUT F 

FRINT F 

IF F=1 THEN GOTO 310 
IF F=2 THEN GOTO 250 
IF F=3 THEN GOTO 5 
GOTO 530 
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10. Readable Scales for 
Computer Plots 


PROGRAM DESCRIPTION 


Every good plotting program should have the ability to automatically calculate 
scales for the X, Y, or Z axes. If the data to be plotted is generated by the com- 
puter, then the user cannot always anticipate the scales needed. The Timex/Sin- 
clair 1000 can, however, automate the task and, moreover, do it intelligently, to 
come up with ‘‘readable’’ scales. The program, given the range of a variable, pro- 
duces a readable linear scale with uniform intervals. Output of the program can 
then be used to drive any plotting routine. A readable linear scale is defined here as 
a scale with interval width the product of an integer power of 10 and 1, 2, or 5, and 
scale values that are integer multiples of the interval width. 
The above definition permits scale values such as: 


—105, 0.0, 0.3, 14)... 
1.24, 1.26, L28,... 
100, 200, 300, . . «5 etc., 


but would prohibit the following examples: 


—1,4,9,... 
LA. 125,205, es 
О, 4, 8,. . ., etc. 


Given a minimum MIN, a maximum MAX of the array to be plotted, and an 
approximate number of desired scale intervals N, the program computes four 
parameters, a new minimum MINP, a new maximum МАХР, an interval width 
DIST, and a number of scale intervals NP which will best fit the plot and approx- 
imate N. The computed parameters, besides producing a readable scale, will 
satisfy the following inequalities: 


(MIN — DIST) < MINP<=MIN, 
МАХ<-МАХР < (MAX + DIST) and 
N/V2.5 < NP < (NV2.5 + 2). 


A unique feature of the program is that it introduces a narrow gate + J (line 200) 
around the MIN and MAX values to avoid an unnecessarily large range between 
MINP and MAXP due to computer roundoff. 
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INSTRUCTIONS 


Initialize the program by pressing the RUN key, then follow prompts for the 
minimum and the maximum of your array and for the approximate desired 
number of scale intervals. The program then returns a new minimum, a new max- 
imum, the interval size, and the number of intervals. 


EXAMPLES 


Problem Definition 


Find the scale values for an array with MIN, MAX, and М, respectively equal to: 
1. —3.1, 11.1, 5 
2. 345, 510, 20 
3. —1100, — 520, 50 
4. 17.5, 17.9, 10 


Sample Run 


SCALING PROGRAM 


ENTER МІМІМИМ?-3.1 

ENTER MAXIMUM?T11.1 

ENTER APPR. М0. OF INTERVALSTS 
NEW MIN= -4 

NEW MAX= 12 

NEW INTERVAL= 2 

NO. INTERVALS= 8 


ENTER MINIMUM? 345 

ENTER МАХІМИМ?510 

ENTER APPR. NO. OF INTERVALST20 
NEW MIN= 340 

NEW MAX= 210 

NEW INTERVAL= 10 

NO. INTERVALS= 17 


ENTER МІМІМИМ?-1100 
ENTER MAXIMUM?- 320 
ENTER APPR. МО, OF INTERVALSTSO 


NEW MIN= -1100 
NEW МАХ= -920 
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NEW INTERVAL= 10 
NO. INTERVALS= 58 


ENTER МІМІМОМ?17, 9 
ENTER MAXIMUMT?T17.5 
ERROR: MIN>=HAX 


ENTER MINIMNUMNTI7.S 
ENTER MAXIMUM?17.9 
ENTER APPR. М0. OF INTERUALS?10 


NEW MIN= 17.5 
NEW MAX= 17.9 
МЕШ INTERVAL= .05 


NO. INTERVALS= 8 


Discussion of Results 


The four examples show that the program produces ‘‘readable’’ scales for a wide 
range of inputs. The first example would result in scale values of —4, —2, 0, 2, 
..., 12. The second example would yield 340, 350, ..., 510. The third example 
would result in — 1100, — 1090, ..., — 520. The fourth example would result in 
17.5, 17.55, ..., 17.9. The fourth example also shows the error message when the 
minimum and maximum were entered by mistake in reverse order. 


PROGRAMMING REMARKS 


Error checking occurs in line 100. Determination of scales is done in lines 240-270. 
The PRINT statement in line 340 also performs computation to save code. 


PROGRAM LISTING 


20 CLS 

30 PRINT "SCALING PROGRAM" 
40 FRINT 

45 PRINT "ENTER MINIMUM?" 5 
50 INPUT X 

60 PRINT X 

70 PRINT "ENTER MAXIMUM?" ; 
80 INPUT Y 

90 PRINT Y 

100 IF Ү>Х THEN GOTO 150 
110 PRINT "ERROR: МІМ>-МАХ" 
120 GOTO 40 


130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
3350 
240 
290 
260 
270 
280 
290 
300 
310 
320 
330 
540 


345 
360 


READABLE SCALES FOR COMPUTER PLOTS 


PRINT "ENTER APPR. NO. OF INTERVALS?" $ 
INPUT N 

PRINT N 

IF М>0 AND N-INT ABS М THEN GOTO 190 
PRINT "ERROR" 

GOTO 130 

LET L=( Y—-X)/N 

LET J=2/1ES 

LET E=INT (LN HU/LN 10) 

LET Ғ-П/10ЖЖЕ 

LET Vz10 

IF Ғ<5ӘК 50 THEN LET V-5 

IF F«SQR 10 THEN LET V=2 

IF F“SQK 2 THEN LET V=1 

LET C-UX10XXE 

LET G=INT (X/C) 

IF ABS (G+1-X/C)<J THEN LET G=G+1 
LET А-СжС 

LET H=INT (Y/UCOo41 

IF ARS (Y/C+1-H)<J THEN LET H=H-1 
LET B=CxXH 

PRINT "NEW МІМ-"»,А,"МЕШ MAX=",Es 
"NEW INTERVAL=" , 

FRINT С»,"М0, INTERVALS=";sH-G 
GOTO 40 
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Originally published as **Pocket Computer Scales Computer Generated Plots." 
Reprinted from Electronics, December 15, 1981, copyright € McGraw-Hill, Inc. 
1981. АП rights reserved. 








PROBABILITY AND 
STATISTICS 


11. Combinations, 
Permutations, and Factorials 


PROGRAM DESCRIPTION 


This program evaluates the three most important functions in statistics and com- 
binatorics. 

Combinations. The program computes the number of possible combinations of 
N different objects taken from a set of M objects. Changing the sequence of the 
objects in a combination does not result in a new combination. 

The resulting function C(M,N) is also called the binomial coefficient and can be 
represented as follows: 


Formula 11-1 








_м__ м _ 
аи IN ~ (M—N)! NI ` 
„кез Des. у Mp3 4, 
"4 S K — Т a Wa 





The first expression, which is the definition of the binomial coefficient C(M,N), 
may lead to large errors by dividing a large number by another large number. If 
any of the factorials in the nominator or denominator exceeds 33! then the first ex- 
pression could not be evaluated by the Timex/Sinclair 1000 at all. On the other 
hand the second expression avoids the overflow problem, hence is the one used 
here to evaluate C(M,N). 

Permutations. The program computes the number of permutations possible 
with M objects taken N at a time. Changing the order of the objects in a permuta- 
tion results in a different permutation. 

The resulting function P(M,N) is computed as follows: 
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COMBINATIONS, PERMUTATIONS AND FACTORIALS 45 


Formula 11-2 


P (M,N) = M (M = 1)....(M = N + 1) 


Factorials: A factorial F(M) equals to the number of ways in which M objects 
can be rearranged. The function is defined and calculated as follows: 


Formula 11-3 


ҚМ) = M (M - 1)....1 


EXAMPLES 


Problem Definition 


1. In how many ways can 10 people be arranged in a row? Evaluate F(10). 

2. How many groups of 3 people can be selected out of a group of 20 when the 
order in each group of 3 matters? Evaluate P(20,3). 

3. How many groups of 3 can be selected out of 20 people, when the order does 
not matter? Evaluate C(20,3). 

4. Same as example 3 but select a group of 17 people out of 20. Evaluate 
C(20,17). 


Sample Run 


FACTORIALS isPERMUTATIONS 2 
COMBINATIONS (BIN. COEFF.) 371 
FACTORIAL ОҒ?10 

FACT( 102-5628800 


FACTORIALS 1»PERMUTATIONS 2 
COMBINATIONS (BIN. COEFF.) 3572 
M-720 

№? 3 

PERM( 20» 3 )=6840 
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FACTORIALS 1»FERMUTATIONS 2 
COMBINATIONS (BIN. COEFF.) 373 
М=?20 

N-73 

СОМЕ 2053271140 


FACTORIALS 1»sPERHUTATIONS 2 
COMBINATIONS (EIN, COEFF.) 373 
М=?20 

М-?17 

COMEX 20,17 )=1140 


РАСТОКІА 5 isPERNUTATIOGNS 2 
COMEINATIONS (RIN, COEFF.) 37 


Discussion of Results 


Example one produces a very large number, examples three and four produce the 
same result. The reason for the equality is that when one chooses a group of 17, 
then the remaining group is three— which is the same as example number three. 


PROGRAMMING REMARKS 


Check of input errors is performed in lines 60 and 130. The same PRINT state- 
ment in line 450 is used for both permutations and combinations by assigning a 
proper name to variable C$. The program makes use of the equality C(M,N) = 


С(М,М — М) and selects in line 250 the fastest computation path. 


PROGRAM LISTING 


20 
50 
35 
37 
40 
30 
60 
70 
100 
110 
120 


CLS 

PRINT 

PRINT "FACTORIALS 1» РЕКМОТАТІОМ5 2" 
PRINT "COMBINATIONS (BIN. COEFF.) 37"; 
INPUT F 

PRINT F 

IF Ғ<>1 AND Ғ<>2 AND F<>3 THEN COTO 30 
GOTO 100ЖҒ 

PRINT “FACTORIAL ПЕТ"; 

INPUT M 

PRINT М 


130 


140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
500 
510 
520 
530 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
479 
480 
490 
495 
200 
210 
313 
520 


530 
540 


COMBINATIONS, PERMUTATIONS AND FACTORIALS 


IF М==0 OR M<>INT ARS M OR M>33 


THEN GOTO 100 

LET В=1 

FOR I-1 TO M 

LET R-BXI 

NEXT I 

PRINT "FACT" МУ" )y=" yE 
GOTO 30 

LET С%-"РЕКМ(" 

GOSUR 480 

IF М>-М THEN GOTO 250 
LET R=0 

GOTO 450 

LET B=1 

FOR I=M-N+1 TO M 

LET Е-ЕЖІ 

NEXT I 

GOTO 450 

LET C$-"CDME(" 

GOSUR 480 

IF М>-М THEN GOTO 350 
LET Е-0 

GOTO 450 

IF М<>М THEN БОТО 380 
LET В=1 

GOTO 450 

LET YzN 

IF М<2ЖМ THEN LET Ү-М-М 
LET К=М+1 

LET R-1 

FOR 1=1 TO Y 

LET Е-ВЖ(К/1-1) 

NEXT I 

PRINT CEMS" SNS” )y=" yE 
GOTO 30 

PRINT "M=?" y 

INFUT M 

FRINT M 

PRINT "М-?"; 

INFUT N 

PRINT N 


IF M=INT ABS M ANI! N=INT ABS N AND МЖҰМ<>0 


THEN RETURN 
PRINT "ERROR" 
GOTO 480 
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12. Error Function ЕВЕ and Its 
Complement ERFC 


PROGRAM DESCRIPTION 


Attributes of many objects in nature, for example the height in inches of the male 
population of a certain city or the population IQs, follow closely the Normal or 
Gaussian distribution—the familiar bell curve. A particular Normal distribution 
is fully described by two parameters, the mean and the standard deviation. The 
Error Function (ERF) is equal to the cumulative Normal distribution and the 
Complementary Error Function (ERFC) equals to 1 — ЕКЕ according to the fol- 
lowing definitions: 


Formula 12-1 


X зашо 
1 f. 1 idt 


ЕП (jo, X) ----- 
ovin —= 


Erfc = 1 — Erf 


The above formulas cannot be evaluated in closed form. The method used in 
this program to compute ERF and ERFC is iterative and leads to results that are 
correct within the full display capability of the Timex/Sinclair 1000 for both 
functions. By comparison, the standard method used to compute ERF and ERFC 
is based on polynomial expansion of a fifth or sixth degree and is precise to only 
3-4 decimal places. ERFC computed via polynomial expansion completely loses 
accuracy at the tails of the distribution. 


INSTRUCTIONS 


Pressing the RUN key issues prompts for the distribution mean, standard devia- 
tion, and the upper distribution limit at which ERF and ERFC should be 
evaluated. 
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EXAMPLES 
Problem Definition 


1. Given a population with a mean IQ of 100 and a standard deviation of 15, 
what fraction of the population will have IQs larger than 135? 

2. The average weight of a bag of fertilizer is 45-165. The manufacturing process 
introduces a standard deviation of 1.5-Ibs. What ratio of bags will be heavier than 
49.5-]bs? 

3. You have probably heard about the 2.33 sigma limit (sigma is the standard 
deviation). What percentage of population does it include? 


Sample Run 


ERROR FUNCTION FROGRAM 

DISTR. MEAN=7100 

ST. DEVIATION-715 

UPPER DISTR. LIMIT-7135 
ERF=0.99018467 ERFC=.0098153287 


DISTR. МЕАМ-?45 

ST. DEUIATION-71.5 

UPPER DISTR. LIMIT-2749.5 
ERF=0.9986501 ERFC=.001349898 


DISTR. MEAN-7100 

ST. DEUVIATION-710 

UPPER DISTR. LIMIT=7123.3 
ERF=0-99009693 ERFC=.0099030756 


Discussion of Results 


Only 0.98% of the population will have IQs larger than 135. The second example 
shows that the danger of heavy bags is rather small, only 0.13%. The third exam- 
ple shows that 2.33 standard deviations include 99% of the population. 


PROGRAMMING REMARKS 


Depending on the value of the argument, the functions are evaluated as continu- 
ing fractions with a fixed number of elements or are evaluated iteratively until they 
converge. 
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PROGRAM LISTING 


20 

30 

35 

40 

30 

60 

70 

80 

90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
560 
370 
380 
590 
400 
410 


CLS 

FRINT 

FRINT "ERROR FUNCTION РКОСКАМ" 
РКІМТ "DISTR. MEAN-7"; 
INPUT R 

PRINT R 

PRINT "ST. DEUIATION-27?"; 
INFUT C 

FRINT C 

IF C>=0 THEN GOTO 130 
PRINT "ERROR" 

GOTO 70 

PRINT "UPFER DISTR. LIMIT-2?"; 
INFUT I! 

FRINT Il 

LET Х-(П-БО/С 

LET X=X/SQK 2 

LET Y-z1/50KR FIXEXF (-XXX) 
IF X<SQR 2 THEN GOTO 280 
LET A=14/X 

FOR I-27 TQ 1 STEP -i 
LET A=I/2/( Х+А ) 

NEXT I 

LET Z=Y/( X+ñ ) 

LET 7=7/2 

LET М-1-7 

GOTO 390 

LET 1=1 

LET Tz2XXXY 

LET N=T 

LET Tz2XXXXXT/€€2XIt41) 
LET M=N+T 

IF M=N THEN GOTO 370 

LET М-М 

LET I=I+1 

GOTO 310 

LET NzN/2t.5 

LET 7=1-М 

FRINT "ERFz"sN»"ERFC-";Z 
FRINT 

GOTO 40 


13. Binomial and 
Hypergeometric Distributions 


PROGRAM DESCRIPTION 


These two distributions occur in many problems involving sampling theory, qual- 
ity control, and games. The binomial distribution applies to the case of a large uni- 
verse, while the hypergeometric distribution applies to cases of a small universe, 
where the samples are taken with no replacement. For example, the probability of 
finding X defective items in a sample from a large universe with a known propor- 
tion of defective items, would follow the binomial distribution. Similarly, the 
probability of finding exactly two red balls in a sample of five balls taken all at 
once from a universe of ten balls, three of which are red, would follow the hyper- 
geometric distribution. Both distributions are very useful for estimating real life 
events. The following formulas apply for the two distributions: 


A. Binomial distribution. 
Given a large universe with a proportion of P special items (defective, marked, 


etc.) and a sample size of N, the probability of finding X special items in that sam- 
ple is as follows: 


Formula 13-1 


Prob (N,P,X) = = р“ (1 — py 


B. Hypergeometric distribution. 
Given a small universe of size M with K defective items, the probability of finding 
X defective items in a sample of size N is as follows: 


Formula 13-2 
K\ /M — K 
X] IN — X 
Prob (M,N,K,X) = м 
М 
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INSTRUCTIONS 


The first prompt after pressing the RUN key gives the choice of the binomial (1) or 
hypergeometric distribution (2). If the binomial distribution is selected, there 
follows a prompt for the sample size N and for the proportion of special items P in 
the universe. In this and following examples, when the prompt calls for more than 
one input value, key in each value separately followed by ‘‘Enter.’’ The next 
prompt is for the number of special items X in the sample for which the probability 
should be computed. After displaying PROB (X) and cumulative probability, the 
next prompt gives a choice between a new X or a new run. 

Similarly, for the hypergeometric distribution, the prompting sequence is the 
same, except that input of three parameters M, K, and N (size of the universe, num- 
ber of special items and the sample size) is requested to describe the distribution. 


EXAMPLES 


Problem Definition 


1. Binomial distribution. It is known that 4% of a manufacturer’s brand are 
defective. What is the probability of finding at most one defective unit in a sample 
of five? 

2. Hypergeometric distribution. In a box of 20 chocolates, five are dark, 15 are 
light. What is the probability of finding one or more dark chocolates in a handful 
of three? 


Sample Run 


BINOMIAL 1» НҮРЕКСЕОМ, 271 


ENTER МэР?5 04 

Х=? 1: 

РКОВ( 1 )=0., 169862931 

SUM OF Р‹Х )= 0.169869 31 
NEU X 1»NEW RUN 271 

X770; 

РКОВ‹ O )=0 ,8153727 

SUM OF FC X)= 0.98524201 


NEW X 1»NEU RUN 272 


BINOMIAL 1»ЭНҮРЕКСЕОМ. 272 


ENTER МэК»М 20 
J 3 
Х-?1 


FROR 20»5»3,19- 0.46052632 
SUM OFC Х)-0.46052652 
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NEU X 1»ЭМЕМ DISTR. 271 
Х=?2 

РКОВ( 20 5» 3,02 )= 0.13157895 
SUM OFC X)=0.59210526 

NEW X 1»МЕЫ DISTR. 271 
Х-?3 

PROR( 20,5»3»30- .008771929S 
SUM OFC Х)-0.60087719 

NEW X lsNEW DISTR. 271 
Х-?0 

РКОБ( 205» 3,Оо- 0.39912281 
SUM OFC X)=1 

NEW X 1,МЕШ DISTR. 2? 


Discussion of Results 


In the first example, the probability of finding one defective unit in a sample is 
0.169. ., and of finding no defects is 0.815... The probability of finding 0 or 1 de- 
fective items in a sample of 5 is therefore P(0) + P(1) = 0.985... 

In the second example the probability of finding one or more dark chocolates in 
a sample of size three is the sum of probabilities of finding one, two, or three, and 
is equal to 0.6. The same result can more easily be obtained by computing the 
probability of no dark chocolates and subtracting it from 1.0, as shown in the last 
run of the example. 


PROGRAMMING REMARKS 


Extensive checking of input values is performed in lines 90, 150, 190, 280, 360, 
and 400. Binomial coefficients are computed in the subroutine between lines 520 
and 620 in the same way as in Program 11. Note use of PRINT statements in lines 
230 and 460, which includes reprinting of input parameters for easier interpreta- 
tion of results. 


PROGRAM LISTING 


20 CLS 

40 FRINT 

41 PRINT “BINOMIAL 1ЭНҮРЕКСЕОМ. 27"; 
45 LET Z-0 

50 INPUT W 

60 PRINT W 

70 IF W-1 THEN GOTO 100 

80 IF W=2 THEN GOTO 290 
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90 
100 
110 
120 
130 
140 
150 


160 
170 
180 
170 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
500 
310 
320 
330 
340 
550 
560 
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570 
380 
590 
400 
410 
420 
422 
424 
426 
428 
430 
432 
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GOTO 40 

PRINT "ENTER МР?" 

INPUT М 

PRINT Nes 

INPUT F 

PRINT Р 

IF N<>INT М OR М<-0 OR Р<=0 OR P>=1 
THEN GOTO 100 

FRINT "X-?"; 

INPUT X 

FRINT X 

IF X<>INT X OR Х<0 OR XN THEN GOTO 160 
GOSUE 520 

LET Т=ЕЖРЖЖХЖ( 1-P OKC N-X ) 

LET Z=Z+T 

PRINT "PROB("$X$")="$T9"SUM OF РХ )=",Z 
PRINT "NEU X i»NEW RUN 27"; 

INFUT W 

PRINT W 

IF Ш=1 THEN GOTO 160 

GOTO 40 

PRINT "ENTER MsKyN"s 

INPUT E 

PRINT Е, 

INPUT G 

PRINT б» 

INPUT F 

PRINT F 

IF E<>INT E OR G<>INT G OR F<>INT F 
THEN GOTO 290 

IF E<=0 OR Е<=0 OR б<=0 OR F>E OR G>E 
THEN GOTO 290 

PRINT "X=?"5 

INFUT D 

PRINT D 

IF П<>ІМТ D OR D<O OR DG THEN GOTO 370 
LET N=G 

LET X-D 

GOSUE 520 

LET J-R 

LET N=E-G 

LET Х-Ғ-П 

GOSUR 520 

LET J-JXh 


434 
436 
438 
440 
450 
460 
465 
470 
480 
490 
500 
510 
520 
530 
540 
550 
550 
570 
580 
590 
600 
610 
620 
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LET N-E 

LET X-F 

GOSUR 520 

LET J-J/R 

LET Z=Z+.J 

PRINT "PROR("sEs"»"5s6G5s","5 
PRINT Fs"s»"sIs5")z",J,"SUM OFCX0z" 97 
PRINT "МЕМ X 1»NEU DISTR. 27"; 
INPUT W 

PRINT W 

IF W=1 THEN GOTO 370 

GOTO 40 

LET YzX 

IF X>€N-X) THEN LET YzN-X 
IF Ү<>0 THEN GOTO 570 

LET R-1 

RETURN 

LET B-1 

LET L=N+1 

FOR I-1 TO Y 

LET B=RK(L/I-1) 

NEXT I 

RETURN 
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14. Uniform and Normal 
(Gaussian) Pseudo-Random 
Deviates 


PROGRAM DESCRIPTION 


This program produces a sequence of either uniform or normally distributed 
(Gaussian) random deviates. Such sequences are frequently used for modeling 
with Monte Carlo method where a direct computation is not feasible. 

An example of an application that is not amenable to exact calculation is that of 
telephone customer satisfaction. The satisfaction of a group of telephone sub- 
scribers depends on the probability distribution of talker loudness, distribution of 
noise, and distribution of transmission losses on a telephone connection. Though 
each use of these distributions can be estimated as Gaussian with a mean and stan- 
dard deviation known from surveys, finding the combined distribution is com- 
putationally difficult. However, the combined effect of all these factors can be 
easily simulated by repeatedly generating pseudo-random deviates from each of 
these distributions and multiplying the individual probabilities. 

This program would normally function as a subroutine in a user-provided ap- 
plication. The formulas used to generate the uniform and normal pseudo-random 
deviates are as follows: 


A. Uniform Distribution 


Formula 14-1 


X, = Fract (997 X ,, + T) 


where RND is the Timex/Sinclair 1000 BASIC pseudo-random function, and X 
and Y are thelower and upper bounds of the random variable. The above formula 
generates uniform distribution values between the lower and upper bounds X and Y. 


B. Normal (Gaussian) Distribution 


Formula 14-2 
X, = (= 2 I U)” sin (27 Ú.) 
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where U, and U, are two sequentially generated uniformly distributed pseudo-ran- 
dom deviates between 0 and 1. The resulting normally distributed pseudo-random 
deviates have a mean of 0 and a standard deviation of 1. The values are then trans- 
formed by the program into a distribution with an arbitrary mean M and standard 
deviation S by multiplying each random deviate with S and then adding M. 


INSTRUCTIONS 


Press the RUN key to start the program, then follow prompts for uniform (1) or 
Gaussian (2) deviates. The next prompt is for the starter which initializes the 
Timex/ Sinclair 1000 RND function. The starter value should be between 0 and 
65535. Entering 0 will always start a different random sequence; otherwise the 
sequence is predetermined. Then follow prompts for the lower and upper bounds 
for the uniform distribution, or for the mean and standard deviation of the Gaus- 
sian distribution, respectively. 


EXAMPLES 


Problem Definition 


1. Generate a number of uniformly distributed deviates between 10 and 20. 

2. Generate Gaussian deviates with mean of 100 and standard deviation of 10. 
A possible interpretation of the resulting values of the Gaussian deviates could be 
a random reading of IQs. 


Sample Run 


UNIFORM 1, GAUSSIAN 271 
STARTER=?123 

LOWER» UPFER BOUNII?T10; 20 
RANDOM ПЕЧІАТЕ=11, 418915 


MORE VALUES 1» NEW RUN 271 
RANDOM DEVIATE=16, 428375 


MORE VALUES 1» NEW RUN 271 
RANDOM DEUIATE-12.132111 


MORE VALUES 1» NEW RUN 271 
RANDOM DEVIATE=19. 917297 


MORE VALUES 1, NEW RUN 271 
RANDOM DEVIATE=13.797302 
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MORE VALUES 1» NEW RUN 271 
RANDOM ПЕУТАТЕ-14.804688 


MORE VALUES 1, NEW RUN 272 


UNIFORM 1» GAUSSIAN 272 
STARTER=?S5S5 

MEAN=T100 

ST. ПЕЧ, =?10 

RANDOM DEVIATE=90. 646835 


MORE VALUES 1» NEW RUN 271 
RANDOM DEVIATE=114.96507 


MORE VALUES 1» NEW RUN 271 
RANDOM DEVIATE=108,. 74587 


MORE VALUES 1» NEW RUN 271 
RANDOM DEVIATE-116.17717 


MORE VALUES 1» NEW RUN 271 
RANDOM DEVIATE=100,. 46652 


MORE VALUES і» NEW RUN 271 
RANDOM DEVIATE-114.597181 


Discussion of Results 


The uniform deviates are between 10 and 20, the Gaussian deviates follow the nor- 


mal distribution. 


PROGRAMMING REMARKS 


Error checking of all inputs is performed in lines 130, 220, and 280. The program 
shows the use of the BASIC RAND and RND functions. If the starter value, or 
argument, of RAND is selected as 0 then a new random sequence will be started 
every time the program is run. Otherwise the pseudo-random values can be 
repeated, if the same starter is chosen, in successive runs. 


PROGRAM LISTING 


20 CLS 
25 FRINT 


UNIFORM AND NORMAL (GAUSSIAN) PSEUDO-RANDOM DEVIATES 


PRINT 
INPUT 
PRINT 
PRINT 
INPUT 
PRINT 


"UNIFORM 1, 
В 

В 
"ОТАКТЕК=?" ; 
5 

5 


GAUSSIAN 27"; 


RAND (5) 
IF B-1 THEN COTO 80 
IF B-2 THEN GOTO 230 
GOTO 30 


PRINT 
INFUT 
FRINT 
INPUT 
PRINT 


"LOWER » 
X 
Xs 
Y 


"өзүү 


UPPER BOUND?" ; 


IF Ү>Х THEN GOTO 160 


PRINT 


"MIN>=MAX" 


GOTO 80 


PRINT 
FRINT 
PRINT 
INPUT 
FRINT 


"RANDOM BEVIATE-Z"SRNIUKCY-XO4X 


"MORE VALUES 1» 
Z 
Z 


NEW RUN 27"; 


IF 2-1 THEN СОТО 160 


IF Z=2 THEN GOTO 25 

GOTO 170 

FRINT “MEAN=?" ; 

INFUT C 

PRINT С 

PRINT "ST. ПЕМ, =?" ; 

INFUT I 

PRINT It 

IF П>-0 THEN GOTO 310 

FRINT "ST. ПЕЧ, < 0" 

GOTO 255 

FRINT "RANDOM DEVIATE-"; 

PRINT SQE (-2XLN RNIOXSIN (2XFPIXKNID)OXIMC 
FRINT 

FRINT "MORE VALUES 1» NEW RUN 27"; 
INFUT Z 


PRINT 


7 


IF 2-1 THEN GOTO 310 
IF Z=2 THEN GOTO 25 
GOTO 320 
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MATHEMATICS 


15. Fourier Coefficients 
of Periodic Functions 


PROGRAM DESCRIPTION 


A continuous periodic function y = f(x) with a period T defined between — T/2 
and + T/2 can be represented as the sum of sine and cosine terms according to the 
following classical formulas attributed to the French mathematician Fourier: 


Formula 15-1 





i A, “ Ao x T 
f(x) = 5 - (A, cos2n т x + B, sin 2n = х) = 25,25, Cn cos (2n Т х+Ө) 
/? > 
A, = = J Қа) cos =F dx C = \/ A24 В? 
=T 
) T/2 А Вп 
2n T X 

= Ө = Arctan(— —— 

B. T f(x) sin T dx ( rm ) 


This program evaluates the above expressions by performing numerical integra- 
tion at specified sampling intervals. If it is known in advance that the function is 
even, f(x) = f(— x), or that it is odd, f(x) = —f(—x), then the program can be 
speeded up because only A or B coefficients have to be evaluated and the integra- 
tion can proceed over only one half of the period T. The program takes advantage 
of this information. 

Though the program only displays the successive Fourier coefficients, these co- 
efficients could easily be stored at the time of computation in a dimensioned 
variable (memory size permitting) for future use. 


60 


FOURIER COEFFICIENTS OF PERIODIC FUNCTIONS 61 


INSTRUCTIONS 


Before running the program, write a subroutine defining the function to be 
evaluated starting on line 1000 and ending with a RETURN (see program listing 
and example). Pressing the RUN key initializes the program and issues a prompt 
for the type of function: (1) for an even function, (2) for an odd function, and (3) 
if the type of function is not known. The subsequent prompts request the size of 
the period T and the number of sampling intervals (two or more). Selection of a 
large number of intervals results in greater accuracy but increased computation 
time. For an even or odd function the requested number of intervals applies to 
only one half of the period T, thus making the intervals smaller and the accuracy 
of the computation better. When a Fourier coefficient is found the computer will 
flash its value. Successive coefficients are computed and displayed until the screen 
is full (error 5) or the BREAK key is pressed. 


EXAMPLE 


Problem Definition 


The function specified in lines 1000-1020 is shown in Figure 15-1. Its period is be- 
tween — 3.5 and +3.5. The program computes successive Fourier coefficients AO, 
Al, A2.... for this even function of x. The approximated curve is shown dashed 
in the same figure. The computation is to be performed with 10 and 100 intervals. 


Figure 15-1 
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Sample Run 


EVEN FUNCT. 1/000 2; 
NEITHER 371 
РЕКТОП Т-?/ 

NO. INTERUALS-710 
A0=4,025 

ñl=3, 0177238 
A2Z=1.1368327 
A3=0.1117739 
А4-.070788988 
AS=0.14849242 
Аб-.058580645 
А/-»01478557 
А8-.027058956 
А9-.011989939 
А10-0 
А11--.011989941 
А12---.027038986 
615-—.014785568 
Al4=—.058580644 
AlS=-0.14849242 
А16--.070788986 
А1/--0.1117739 


Discussion of Results 


EVEN FUNCT. 1>0ШІ 2; 
NEITHER 371 
РЕКІПІ T=?7 

NO. INTERVALS=?100 
А0-4.00015 
A1=3.0345888 
А2-1.1793983 
Аӛ-0.10390619 
А4- 058466219 
А5-0.18842451 
Аб-.08414665 
А/-»0001493158 
АВ8-.047250562 
А9-.0579352911 
А10-.0094089601 
А11-.0077683475 
AL2=.0324 3459 
А13-.017811636 
А14-.00014726251 
А15-.013311727 
A16=.018110545 


The Fourier coefficients were independently verified with a high precision algo- 
rithm on a different computer; it was found that with 100 intervals between four 
and five decimal places on the Timex/Sinclair 1000 are correct. Note the Gibbs 
phenomenon at sharp corners, X = —3.5,0, +3.5; the Fourier expansion does 
not converge at those points to its true value. For example, at X = 0 the series 
seems to converge toward 6.85 rather than 7.0 as expected. At other points the 


convergence is excellent. 


PROGRAMMING REMARKS 


Checks on input errors are performed in lines 60, 100, and 140. The program per- 
forms trapezoidal integration of the function starting in line 1000. 
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PROGRAM LISTING 


30 

40 

930 

60 

70 

80 

90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
195 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
500 
510 
320 
530 
340 
350 
455 
560 
370 
380 
390 
400 
410 


PRINT “EVEN FUNCT. 1,ОПП 2,МЕТТНЕК 37"; 
INPUT А 

PRINT А 

IF А<1 OR A>3 OR A<>INT A THEN COTO 30 
PRINT "РЕКІОП Т=?"; 
INPUT В 

РКІМТ В 

IF В<=0 THEN GOTO 70 
PRINT "NO. INTERVALS=?"s5 
INFUT C 

PRINT C 

IF С<2 OR C<>INT C THEN GOTO 110 
IF A=3 THEN GOTO 500 

LET [i=B/2/C 

IF A=2 THEN GOTO 340 

LET J=-1 

LET J=J+1 

PAUSE 20 

LET X=0 

GOSUR 1000 

LET EzY 

LET F=0 

LET G-2XJXFI/RB 

FOR 1=1 TO C 

LET Х-ПЖІ 

GOSUR 1000 

LET Ғ-ҒЖК ЕҒҮ)/2ЖПЖСО5 (GkC X-II2)0) 
LET E-Y 

NEXT I 

LET H-4XF/B 

PRINT "A" Ју "=" Н 

GOTO 190 

LET J=0 

LET J=J+1 

PAUSE 20 

LET X=0 

GOSUR 1000 

LET E=Y 

LET K=0 

LET G=2x.JKXFI/E 

FOR І-1 TO C 
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420 LET X=LIXKI 

430 GDOSUR 1000 

440 LET K-N*CETYO/2XIXKSIN (CKO Х-П/2)) 
450 LET E=Y 

440 NEXT I 

470 LET L=4xK/B 

480 PRINT "B" yJy"=" yL 

490 GOTO 350 

900 LET D=B/C 

310 LET J--1 

520 LET J=J+1 

525 PAUSE 20 

350 LET X--B/2 

540 GOSUB 1000 

290 LET ESY 

960 LET K=0 

370 LET F=0 

580 LET б-2ЖУЖРІ/Е 

390 FOR I-1 TO С 

600 LET X-DXI-B/2 

610 GOSUR 1000 

620 LET М-бЖ( Х-П/2) 

630 LET F=F+( E+Y 0/2XIXXCOS М 
640 IF J-O THEN COTO 660 
650 LET К=К+СЕ+Ү )/2x0KSIN M 
660 LET E=Y 

670 NEXT I 

680 LET Н-2ЖҒ/В 

690 PRINT "“A"sJs"="3H 

700 IF J-0 THEN GOTO 520 
710 LET Lz2XK/R 

720 PRINT "E"sJ5"z";L 

730 LET N=SQR (HXH4LXL ) 

740 LET Р-АТМ C-L/H) 

750 PRINT "С" ЈУ "=" М 

760 FRINT "ТНЕТА" ЈУ" =" F 
770 GOTO 520 
1000 LET Ү-7-3.5ЖАВ5 X 
1010 IF ABS Х>2 THEN LET Y=0 
1020 RETURN 


16. Runge-Kutta Solutions 
of Differential Equations 


PROGRAM DESCRIPTION 


Runge-Kutta methods for numerical solution of ordinary differential equations 
consist of a group of iterative formulas distinguished by certain important com- 
mon features: 

1. The derivative y’ can be given as a function of x and y. It does not have to be 
provided as an explicit function of x alone. 

2. The method is self-starting. Knowledge of the function is only required at 
one point to compute the succeeding point in the iteration. 

3. The integration interval can be changed as the iteration proceeds. 

4. The formulas usually converge quickly. 

Successively higher order Runge-Kutta formulas and finer integration intervals 
lead to more precise solutions. The program evaluates fourth order Runge-Kutta 
formulas for solving first order differential equations of the type: 


у' = f(x,y) 
The formulas assume that f(x,y) is known. The initial condition X0 and YO, the in- 


tegration interval H along the X-axis, and the final value of X are selected by the 
user. The program evaluates the following formulas: 


Formula 16-1 


K, Hf(X;, Y) 


, = НХ, +Н/2, Y. +K 


K 
К, = HEX, H2, Y, + К //2) 
К 


4 = Hf(X;+H/2, Y; + K/2) 


| 


ү Yee (K, +2K,+2K,+K,) 
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INSTRUCTIONS 


Before running the program write the function y’ = f(x,y) as a BASIC subroutine 
starting on line 1000 and append it to the program as shown in the program listing 
and in the example. Use the letter R for the derivative. Pressing the RUN key starts 
the program and prompts for the interval size Н, the initial values of ХО and YO, 
and for the final value of X at which the solution should be found. The interval H 
should be positive, though X can be smaller or larger than X0. There are no other 
limitations on H. The program flashes the succeeding x values as they are 
evaluated and stops when the computation is completed by displaying Y(X). 


EXAMPLE 


Problem Definition 


Given the differential equation y’ = x + y, with the initial condition X= 0, Y = 0. 
Find Y(1). Assume integration intervals of 0.15 and 0.05. 


Sample Run 


4-TH ORDER RUNGE KUTTA 
INTEGRATION INTERVAL H?0.15 
ENTER XO»YO»XLASTTOs O51 
0.15»0.3»,0.45»0.6»0.75,0.9,1» 
Ү(12-0.71827251 


INTEGRATION INTERVAL H7.05 

ENTER X0OrYOsXLAST?TOs021 
.05»0.1»,0. 15»0.2»0. 25:0. 5»0. 3550 
.4,0.45,0.5,0.55»0.6»0.65»0.7»0. 
753»0.8»0.85»,0.9»>0.95:1,1», 

Ү( 12-0.71828169 


Discussion of Results 


This happens to be a differential equation for which an exact solution can be 
found. The solution is: 


ҮС) =e — 2 = 0.718281... 
The Runge-Kutta solution is 0.718272... for Н = 0.15 and 0.718281 for the inte- 


gration interval of 0.05. The shorter integration interval evidently produces a 
more exact solution. 
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PROGRAMMING REMARKS 


Input error check on the interval size is performed in line 60. Notice the display of 
the final value of Y(X) in line 245. The statement uses XLAST as its argument for 
clarity. As mentioned in the Program Description, the routine starting in line 1000 
describes the particular differential equation to be evaluated. 


PROGRAM LISTING 


20 PRINT "4-TH ORDER RUNGE KUTTA” 
30 FRINT "INTEGRATION INTERVAL Н?" ; 
40 INFUT A 

30 PRINT А 

60 IF A<=0 THEN COTO 30 
70 PRINT "ENTER XO»YO»XLAST?"s 
80 INPUT B 

90 PRINT В; 

100 INPUT C 

110 PRINT ",“%С; 

120 INPUT I! 

130 PRINT ","%П 

135 IF В=П THEN GOTO 70 
140 LET E=INT ABS (П-ВО/А 
150 LET Ғ-Е 

160 LET G=SGN (1-Е )ЖА 

170 LET H=C 

180 FOR I-1 TO E 

190 GOSUR 270 

200 NEXT I 

210 IF Ғ-П THEN GOTO 240 
220 LET G=LI-F 

230 GOSUB 270 

240 FRINT 

245 PRINT "Y(" yI" y=" +H 
250 PRINT 

260 GOTO 30 

270 LET X-F 

280 LET Y=H 

290 GOSUE 1000 

500 LET K=GXR 
310 LET X=F+G/2 
320 LET Y=HtK/2 
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220 
340 
350 
340 
370 
380 
390 
400 
410 
420 
430 
440 
450 
440 

1000 

1010 
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GOSUE 1000 
LET L-CXR 
LET Y-HtL/2 
GÜSUR 1000 
LET M-GXR 
LET X=F+G 
LET Y=H+M 
GOSUR 1000 
LET N-GXR 
LET H=HtK/6+L/3+M/34N/6 
LET F-F4G 
PRINT Fj","j 
PAUSE 20 
RETURN 

LET К=Х+Ү 
RETURN 


17. Numeric Solution of 
Transcendental Equations 


PROGRAM DESCRIPTION 


This program solves numerical equations which in general cannot easily be solved 
explicitly. For example, it will solve for X: 


F(X) = LOG X - X^ + 10=0. 


The program requests an initial estimate of ХО in the neighborhood of the solu- 
tion. Such an estimate can usually be made from a plot of the function. It then 
requests a DELTA increment at which the function will be evaluated (a large 
DELTA will converge more quickly, but may miss the true solution). The program 
then requests EPSILON, a small positive number. When the absolute value of the 
function is less than EPSILON, iteration stops and the solution is printed. If more 
than one solution exists, then the program will find the solution closest to the ini- 
tial estimate. 

The program uses a modified Newton's iteration method to find the zero cros- 
sing of the function, which is the required solution. The program starts by **draw- 
ing” a straight line through the function evaluated at ХО + DELTA/2 and X0 — 
DELTA/2. The crossing of this line with the X axis gives the next estimate of ХО. 
If the estimate of XO stops changing while the absolute value of F(X) is still larger 
than EPSILON, then the DELTA interval is halved by the program and the itera- 
tion continues. It is obvious that this method will only work for functions without 
irregularities between the first estimate of ХО and the true ХО. The function F(X) 
also has to be single-valued and its first derivative cannot beO in this region, other- 
wise the straight horizontal line will not cross the X-axis. To indicate how the 
iteration is progressing (hopefully it is converging!), the program flashes suc- 
cessive approximations of ХО. As long as these are converging you should arrive 
at the approximate solution within a short time. The program then stops and 
displays its best estimate of the solution to the equation. 


INSTRUCTIONS 


As in the previous program, the equation to be solved is appended to the program 
before running it. The subroutine describing the equation should start at line 1000 
and should end with a RETURN as shown in the listing and example. Pressing the 
RUN key starts the program and generates three prompts, for the initial guess of 
X0, DELTA, and EPSILON. The iteration then continues until ABS (F(X)) 
< EPSILON. 
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EXAMPLE 

Problem Definition 

We will try to solve the equation mentioned earlier: 
F(X) = LOG X - X° + 10 = 0. 


Use the initial guess for ХО = 2.0, DELTA = 0.01, and EPSILON = 1E—6, then 
repeat with a larger DELTA of 0.1. 


Sample Run 


TRANSCENIENTAL EQUATION SOLUTION 
ENTER GUESS FOR ХО, ПЕТА XsEPS.? 
2+ .01,1Е-6 

SUCCESSIVE APPROXIMATIONS? 


3+ 6656754 322673533 
342421 374 34 2420552 
342420352 (24303552 
3, 2420352 3, 2420382 
3, 2420382 Ë, 4429057 
3, 242039 3,242049? 


X0=3, 2420392 


TRANSCENDENTAL EQUATION SOLUTION 
ENTER GUESS FOR ХО, ПЕТА ХэЕР5,? 


2,0.1,1Е-6 

SUCCESSIVE APPROXIMATIONS? 
43, 6650052 5.266929 
342417394 3.2416374 
302416374 3+ 2416374 
3, 2419388 £, 2419388 
3, 2419388 3, 2420141 
5.2420141 5.2420141 
302420329 342420329 
3+ 2420377 £, 2420377 
3+2420388 3+ 2420388 
5+ 2420391 


X075, 2420591 
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Problem Discussion 


Both runs, after happily flashing converging values of X, lead to the result of 
Х-3.242... 


PROGRAMMING REMARKS 


The program flow ends with a STOP in line 260 in order not to collide with the 
subroutine describing the equation. For every run you must press the RUN key. 


PROGRAM LISTING 


150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
245 
250 
260 
1000 
1010 


РКІМТ "TRANSCENDENTAL EQUATION SOLUTION" 
PRINT "ENTER GUESS FOR ХО, ПЕТА ХЭЕР5,?" 
INFUT А 

LET E=A 

РКІМТ А; 

INFUT B 

FRINT ",. "#6; 

INFUT С 

PRINT "y" yC 

PRINT “SUCCESSIVE АРЕКОХІМАТІОМ5; " 
LET Х-А-В/2 

GOSUR 1000 

LET D=Y 

LET X=AtB/2 

GOSUR 1000 

LET Z-BXI//CI!- Y )+A-B/2 

LET X=Z 

FRINT Z: 

FAUSE 30 

GOSUR 1000 

IF ABS Ү<=АЕ5 C THEN СОТО 240 
IF Y=E THEN LET RB-Rh/2 

LET А=7 

LET EzY 

GOTO 90 

FRINT 

PRINT 

PRINT "XOz"5Z 

STOF 

LET Y=LN X/LN 10-XXX110 
RETURN 


18. Complex Arithmetic 


PROGRAM DESCRIPTION 


With this program your Timex/Sinclair 1000 computer becomes a full-fledged 
complex arithmetic calculator operating in Reverse Polish Notation (RPN). For 
those of you familiar with Hewlett Packard calculators, which use RPN, this will 
be old hat. If you are unfamiliar with RPN a few explanations are necessary. 

Think of a stack of registers, four deep, named X, Y, Z, and T in your com- 
puter. The bottom of the stack, the X register and the one immediately above it, 
the Y register, are always being displayed. Each register contains the real and the 
imaginery parts of a complex number. Unary (one variable) operations; e.g., 
squaring, taking the square root, computing reciprocals, etc. are performed on 
the X register alone and do not affect the stack. Operations on two variables, such 
as addition, subtraction, multiplication, and division operate on the X and Y reg- 
isters. The result of a two-variable operation always appears in the X register and 
the stack drops, that is, the contents of register Z are moved into register Y, and 
the register T is copied into Z. Register T retains its previous value. Numbers are 
entered into the stack via the ENTER command (1). Entering a complex number 
moves it into the X register replacing the previous contents of that register. Several 
stack manipulation operations are provided: Х<->Ү exchanges the contents of the 
X and Y registers, RAISE STACK destroys the previous contents of register T, 
moves previous contents of register Z into T, Y into Z, and duplicates X into Y. 
The program also provides a memory register M. Two commands are provided to 
copy register X into M (SAVE), and to copy register M in X (RECALL). 


INSTRUCTIONS 
RUN — Start computation, issue prompts for operations 1-12 
1 —ENTER, enters the complex number into X register, stack remains un- 
changed 


2  —RAISE STACK, copies contents of Z in T, Y in Z, and X in Y; destroys 
former contents of T. 

— Х<->Ү, exchanges the X and Y registers 

— SAVE, copies X into M 

— RECALL, copies M into X 

— Adds X and Y, puts result in X, drops stack 

— Subtracts X from Y, puts result in X, drops stack 

— Multiplies X and Y, puts result in X, drops stack 

— Divides Y by X, puts result in X, drops stack 

— Finds reciprocal of X, puts it in X 

11 —Finds square root of X, puts it in X 

12  — Finds square of X, puts it in X 
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EXAMPLE 


Problem Definition 


Evaluate the following expression: 


Formula 18-1 


(2 + 3j) (5 — 4j) 
5 + 4j 


1 


Sample Run 


OPERATION 1-1271 
ENTER К,І?2»3 


REAL IMAGINARY 
Y: 0 0 
Xi 2 3 
OPERATION 1-1273 
REAL IMAGINARY 
Ж 3 
Хх; О 0 


OPERATION 1-1271 
ENTER К,І?5»-4 


REAL IMAGINARY 
Ti 2 3 
Х% 5 -4 
OPERATION 1-1278 

REAL IMAGINARY 
Үз U 0 
X: 22 7 
OFERATION 1-127343 

REAL IMAGINARY 
ys 22 7 


х: 0 0 
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OPERATION 1-1271 
ENTER К,1?5,4 


REAL IMAGINARY 
Y+ 22 7 
Х% 5 4 
ОРЕКАТТОМ 1-1279 

REAL IMAGINARY 
Y: 0 0 
Хі 3,3638337 -1.2926829 
ОРЕКАТІОМ 1-12711 

REAL IMAGINARY 
Y: 0 0 
Хх, 1.8670036 -0.34619187 
OPERATION 1-12712 

REAL IMAGINARY 
Y: Q 0 
Xs 3,3658537 -1.2926829 


Discussion of Results 


The complex number 2 + 3j is entered first. It is then moved from the X into the Y 
register before the second complex number 5 — 4j is entered. X and Y are then 
multiplied and the intermediate result is 22 + 7j. The product is then moved into 
the Y register and the last complex number 5 + 4j is entered. Next a division and a 
square root operation is performed. The result is then 1.86... —0.346...j. In the 
final step the result is squared to demonstrate that the complex square root and 
squaring are inverse operations. Note that a second root also exists which is the 
negative of the first one, namely — 1.86... - 0.346. . .j. 


PROGRAMMING REMARKS 


The program fills the 2K RAM capacity of the Timex/Sinclair 1000. Therefore 
only the X and Y registers are displayed rather than the whole stack. The square 
root function (11) displays only one root. The other root is the negative of the first 
one. Stack raising and dropping is performed in FOR NEXT loops in lines 100-110 
and 1000-1015. The only input error check is performed in line 27. Multiple 
branching is performed in line 40 by multiplying the operation code by 50. 
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PROGRAM LISTING 


10 
15 
18 
20 
25 
27 
40 
30 
55 
60 
65 
67 
70 
72 


75 


74 

75 
100 
105 
110 
115 
150 
165 
170 
175 
180 
200 
205 
210 
250 
255 
260 
500 
305 
310 
350 
555 
360 
400 
403 
410 


DIM X(8) 
PRINT 

PRINT "OPERATION 1-127"; 
ІМРОТ F 

PRINT F 

IF F<>INT F OR Ғ<1 OR Ғ>12 THEN GOTO 18 
GOTO 50ЖЕ 

PRINT "ENTER RoI?" 
INPUT X(1) 

PRINT X(1); 

INPUT X( 2) 

PRINT ","$X(2) 

LET К-Х(1) 

LET J=X( 2) 

PRINT " REAL"y"IMAGINARY","Y1 
"$3XCS0sXCADS"XE. "SRJ 
LET II-RKR-HJKJ 

GOTO 15 

FOR I-8 TO 3 STEF -1 
LET X( I)-XC I-2) 

NEXT I 

GOTO 70 

LET X(10-X(3) 

LET X(2)=X( 4) 

LET X(3)-R 

LET X(4)=J 

БОТО 70 

LET K-R 

LET M-J 

GOTO 70 

LET Х(1 )=К 

LET X(2)-M 

GOTO 70 

LET X(1)0-R4XC 3) 

LET X( 22-J4XC 4) 

GOTO 1000 

LET X(10-XC3)-R 

LET Х(2 )=Х( 4)-J 

GOTO 1000 

LET X( 1 )=RKX( 3)-JKX( 4) 
LET Х(2)=КЖХ( 4)+JXX( 3) 
GOTO 1000 
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LET XC 102C6RXXCSOTJIXXCAOD/D 


з LET Х(2)=( RKXXC 4 )>JkX( 3))/П 


GOTO 1000 

LET XC 10zR/T! 

LET X( 2»--J/TI! 

GOTO 1000 

IF К<>0 ОК J«-0 THEN GOTO 554 

LET X(1)20 

(ЕТ X( 22-0 

GOTO 70 

LET ХХ190-5ӨК ((RESQR CREXRT JXJ 25/2) 


ч LET X€2)=J/2/X01) 


GOTO 70 

LET X(1)-RXR-JKJ 
LET X( 22-2XRXJ 
GOTO 70 

FOR I-3 TO 6 

LET XCI)=X( 142) 


NEXT I 
GOTO 70 


19. Decimal-to-Binary 
Fractions 


PROGRAM DESCRIPTION 


That part of the civilized world that employs the British system of measurements 
still prefers binary fractions such as 5/16 or 7/32 to numbers as easy to remember 
as 0.3125 or 0.21875. To accommodate both the ‘‘decimal’’ and the ‘‘binary’’ 
population, this program will translate a decimal fraction into its binary equiv- 
alent. We will define a **binary"' fraction as a fraction with a denominator that is 
an integer power of 2. All the program user has to do is to indicate the smallest ac- 
ceptable fraction — a binary power of 2; e.g., 256 for 1/256, 64 for 1/64, etc. 

If the translation can be done exactly, then the program will divide the numer- 
ator by the denominator. For example, 0.375 will be translated to 3/8 though the 
smallest fraction requested by the user may be 1/256. If the translation cannot be 
done exactly, the smallest binary denominator will be found within the given con- 
straints, and the percent error will be indicated. 


INSTRUCTIONS 
Pressing the RUN key issues two prompts: ENTER DECIMAL FRACTION? — 
the number to be translated into the binary fraction, and MAX DENOMINA- 


TOR (PWR OF 2)— the number to indicate the largest acceptable denominator. 
The program then displays the binary fraction and the error, if any. 


EXAMPLE 


Problem Definition 
Translate 2.5, 0.12717, and 0.1875 into binary fractions. Do not accept any 


fraction smaller than 1/16 for the first example and 1/256 for the second and third 
examples. 
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Sample Run 


ENTER DECIMAL FRACTION?2.5 
MAX. DENOMINATOR (РЫК OF 22716 


2.5-5/2 


ENTER DECIMAL РКАСТІОМ?О, 12717 
МАХ. DENOMINATOR (РЫК OF 2)7256 
0.12717 APFROXIMATES 33/256 
WITH ERROR ОҒ -1.3652984 PERCENT 


ENTER DECIMAL FRACTION?TO.1875 
MAX. DENOMINATOR (РЫК OF 227256 
0.1875-3/16 


Discussion of Results 


The first and third examples lead to exact binary fractions. The second example 
сап be approximated wtihin the constraints of 1/256 to within 1.35...%. 


PROGRAMMING REMARKS 


The only error check on input value (largest denominator) is performed in lines 
70-90. To account for the Timex/Sinclair 1000 roundoff a ‘‘fudge’’ factor of 
0.0001 is used in lines 90 and 120. To find the smallest denominator the program 
keeps dividing the numerator and the denominator by 2 (lines 140-150), as long as 
both are integers. 


PROGRAM LISTING 


J PRINT 
10 PRINT "ENTER DECIMAL FRACTION?" ; 
20 INPUT А 
30 PRINT А 
40 PRINT "МАХ. DENOMINATOR (РЫК OF 207%; 
30 INPUT B 
60 PRINT В 
70 IF B«2 THEN GOTO 40 
80 LET C-LN B/LN 2 
90 IF ABS (CC- INT C)>.0001 THEN GOTO 40 


100 
110 
120 


130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
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LET C=INT C 

LET П-АЖВ 

IF (ABS D-INT ABS (/ПІ4.00010)>-.0001) 
THEN GOTO 170 

IF П/2<>ІМТ (072) THEN GOTO 230 

LET D=D/2 

LET B-B/2 

GOTO 130 

LET II-INT (I.5) 

LET Е=А-П/Е 

LET Ғ-100ЖЕ/А 

PRINT As" APPROXIMATES "yI" Z" yE 
PRINT "WITH ERROR OF "$F$" PERCENT" 
GOTO 5 

PRINT А" =" унин yB 

GOTO 5 


OPERATIONS 
RESEARCH 


20. Solution of Queuing 
Equations 


PROGRAM DESCRIPTION 








The program presented here solves the classical M/ M/S queuing equations which 
describe many real life situations. 

The M/M/S queuing equations assume exponential distribution of customer 
arrival rates L, exponential distribution of serving times M, S identical servers, 
and FIFO order of service. Exponential distribution implies independence of 
events; the customers will seek service independent of the queue length and the 
servers will operate at a steady rate independent of the load and queue length. 
Such assumptions usually result in safe (pessimistic) estimates of waiting times 
and other queuing parameters. The M/ M/S cases result in relatively simple equa- 
tions (in terms of queuingtheory) and can be solved within the memory limita- 
tions of the Timex/ Sinclair 1000. The program can be used to estimate such pro- 
blems as how many barber chairs to provide at a midtown location, how many 
telephone operators should be working between 2 and 3am, or how many printers 
should be connected to a central processor. In general, application of the program 
will help in balancing the cost of providing service with the inconvenience of long 
waiting lines. The queuing equations used in the program are given below. 


Program limitations: S < 34, U < = 1.0, all inputs positive, all intermediate re- 
sults should be less than 1.0E38. 


Inputs: 
Av. arrival rate r 
Ау. service rate т 
Number of servers s 
Number of items in system n 


Outputs: 
Facility utilization U =r/(ms) 
Prob. of finding no items in system Р(0) 


P(0) = (с/т) /s!(1-U)+ X (ету ДІ)! 
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Prob. of finding n items in the system P(n) 
P(n) = P(0) x (r/m) x (1/n!) for n<s 
P(n) = P(0) x (r/m)" x (1/(s!s"*)) for n> =s 
Prob. of finding all servers busy B 


B= I P(n) = (r/my) x P(O)/(s! x (1-U)) 


Average waiting in queue Tw — B/(sm X (1-U)) 

Average system response time Tr = Tw + Im 

St. dev. of Tr = Tw/B x V(BQ — B) +s? x (1 — Uy) 

Prob. of Tw>T =B X Exp( - smT x (1 — U)) 

Prob. of Tr^ T = Exp( — mT) X (1 + (1 — Exp( — msTk)) X B/sk) for k <>0 where 
k=1—U-I/s 

Prob. of Tr>T = Exp( — mT) х (1+ BmT) for k=0 

Average number of items in queue Q =r Xx Tw 

Average number of items in system N =r X Tr 


The program provides all prompting and checks input data for correct values. 


PROGRAM INSTRUCTIONS 


Pressing the RUN key will issue prompts to enter the arrival rate L, the service rate 
M and the number of servers S. An incorrect input will cause the prompts to re- 
peat. The program then computes and displays the previously mentioned parame- 
ters of the queuing system, namely B, U, P(0), Tw, Tr, Standard deviation of Tr, 
Q, and N. The program will then issue a prompt “ENTER CODE 1, 2 OR 3?" To 
compute P(n) enter 1, P(Tw, Tr) — 2, or to continue with a new run—3. Selecting 
P(n) will elicit a prompt for the number of items in the system n and will compute 
and display P(n). Choosing P(Tw,Tr) will prompt for time T and will compute 
and display P(Tw » T) and P(Tr » T). The prompts then repeat. 


EXAMPLE 


Problem Definition 


A computer working area has five terminals. On average 16 customers per hour 
arrive at random intervals with each job taking on average 15 minutes (four per 
hour). Compute the queuing parameters; also find the reduction in waiting time if 
a sixth terminal is installed. 


Sample Run 


ENTER 1 »,М»5716»4»5 
PCALL BUSY )= 0.55411255 


UTIL= 0.8 
P( 0 )= ‚012987013 
AV.WAIT= 0.13852814 


AV.RESP.= 0. 38852814 
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ST. ПЕМ, (ТК )= 0. 33552053 
AV. ІМ SYS.= 6.2164502 
АУ IN QUEUE= 2.2164502 


ENTER CODE 1,2 ОК 571 
NO. ITEMS IN SYSTEM?1 
PCN=1 )=, 051948052 


ENTER CODE 1,2 OR 371 
NO. ITEMS IN SYSTEM?6 
FC N26 )=. 088658009 


ENTER CODE 1,2 OR 3571 
NO. ITEMS IN SYSTEM?10 
PC N=10 0=.03631432 


ENTER CODE 1,2 OR 372 
TIME-70.2 

PC TW-0,.2)20,24897882 
PC TR>0. 2 )=0,.64851202 


ENTER CODE 1,2 OR 372 
TIME-7T10 

FC TW7-100722, 3540665E-18 
PC ТК>102)-9.8411012Е-17 


Discussion of Results 


ENTER L»M»S5?16»4»6 
PCALL BUSY )= 0.28476085 


UTIL- 0.66666667 
P( 0 )= „016685206 
AV.WAIT- 035595106 
AV.RESP.- 0. 28559511 


ST. DEV.CTR27 0.26482403 
AV. ІМ SYS.- 4.5695217 
AV IN QUEUE- 0.56952169 


ENTER CODE 1»2 ОК 371 
М0. ITEMS IN SYSTEM?1 
P(N=1 )=,066740823 


ENTER CODE 1,2 OR 372 
ТІМЕ-?0.2 

PC TW-0,2)2.057492223 
FCTR-0,2220,51978804 


The first part of the example (five terminals) shows a fairly high utilization of the 
facility, 0.8. The average waiting time in queue is 0.138 hours or slightly over eight 
minutes. Adding a sixth terminal will cause a considerable improvement. The 
waiting time will drop to 0.035 hours or about two minutes. All other wait prob- 
abilities are also correspondingly affected. 


PROGRAMMING REMARKS 


A check of input values is performed in line 80. Notice use of the roundoff fudge 
factor in line 60. A fairly elaborate PRINT statement is shown in line 900. As you 
enter the program notice the use of the EXP function (Shift ENTER to get in the 
FUNCTION mode, followed by pressing the X key). 
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PROGRAM LISTING 


м 
10 
20 
30 
40 
30 
60 
70 
80 


85 

90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
273 


280 


270 
350 
555 
340 
370 
380 
370 
400 
410 
420 


CLS 

PRINT 
INFUT 
PRINT 
INPUT 
PRINT 
INPUT 
PRINT 


IF L<=0 OR М<-0 OR 5<-0 ОК 5<>ІМТ S OR 8255 
THEN GOTO 10 
LET Р=0 


LET 


X-L/M 


LET UzX/S 

LET V=1-U 

IF У<-0 THEN GOTO 10 
LET Ү=5 

GOSUB 1000 

LET T=Z 

FOR I=0 TO 5-1 
LET Y=I 

GOSUR 1000 

LET PzPtXXXI/Z 
NEXT I 

LET Ү=5 


GOSUB 


1000 


LET Р-1/(РҒХЖЖ5/Т/У) 
(ЕТ B-XXXSXP/T/VU 

LET W-zRB/S/M/V 

LET R-zWt1/M 


PRINT 
PRINT 
PRINT 
kU /E 

PRINT 
FRINT 
PRINT 


"PCALL BUSY )=", В, "ИТІ =", 00у "РО )=", 
Ру "АМ, ШАІТ=" ууу "АИ, RESF. =", К 
"ST.IEU.CTROZ" SRR ( 2KR-BKR+( S*£U )KXK2 ) 
“АУ, 


ІМ бҮб,=",1ЖК,"Л\ IN QUEUE-" LXW 


"ENTER CODE 1,2 OR 357"; 


INPUT K 


PRINT 


IF К<> 


PRINT 
INFUT 
PRINT 


K 

1 THEN COTO 510 

"NO. ITEMS IN SYSTEM?" s 
N 

N 


LET П-РжЖХЖЖМ 
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430 IF М<5 THEN GOTO 460 
440 LET E-D/T/SXX( N-5) 

450 GOTO 490 

460 LET Y=N 

470 GOSUB 1000 

480 LET E-I/Z 

490 PRINT "P(Nz"sN;")z"sE 
400 GOTO 350 

310 IF К<>2 THEN GOTO 5 
920 FRINT "TIME-T7"; 

330 INFUT H 

340 PRINT H 

290 LET F=BKEXP -( 5жМжУЖН) 
360 LET A$="W" 

970 GOSUR 900 

980 LET А%="К" 

370 LET GzU-1/S 

600 LET П=ЕХР -( MXH) 

610 IF ARS G<=.00001 THEN GOTO 640 
620 LET F=DKC 145/8% 1-EXF -CHMXSXGXH)2)0/G) 
630 GOTO 650 

640 LET FzIXC1itTBXMXH) 

650 GOSUE 900 

660 GOTO 350 

900 PRINT "РСТ" АФУ" "НУ" )z" FF 
910 RETURN 

1000 LET 7=1 

1010 FOR Ј=1 TO Y 

1020 LET Z-ZXJ 

1030 NEXT J 

1040 RETURN 


Originally published as ‘‘Pocket Computer Tackles Classical Queuing 
Problems." Reprinted from Electronics, September 22, 1981, copyright © 
McGraw-Hill, Inc. 1981. АП rights reserved. 


21. Reliability Parameters of 
Complex Systems 


PROGRAM DESCRIPTION 


Reliability of a composite system; e.g., a telephone office consisting of incoming 
and outgoing trunks, switching machines, and power supplies or a computation 
center consisting of CPUs, printers, disk drives, and power supplies, can be ex- 
pressed in terms of the Mean Time Before Failure (MTBF), and Mean Time Till 
Repair (MTTR). A high MTBF and alow MTTR can be achieved by using reliable 
components which do not break down frequently and are easy to repair, or by 
providing redundancy. For example, two CPU’s can share the load; if one fails, 
the other takes over. 

The program will find the overall reliability of a system given MTBFs and 
MTTRs of individual components. It is assumed that failures of individual com- 
ponents are independent of one another. The reliability curve then follows the 
center part of the classical **bathtub"' curve, where the high number of failures іп 
the beginning is caused by initial problems and the high number of failures 
towards the end of useful life is caused by component wearout. It is also assumed 
that individual components will be grouped into subsystems of one of three types 
listed below, and that a number of subsystems are connected in series to form an 
overall system. 

A. Series subsystems. Failure of one component will incapacitate the whole 
system. 

B. Parallel subsystems. At least one component has to stay alive for the system 
to operate. 

C. Parallel k-of-n subsystem. Where k components out of n have to stay alive 
for the system to operate. 

The program computes MTBF and MTTR for each subsystem and for the 
overall system. It also finds the availability factor and the probability of x com- 
ponents failing during time t, where x and t can be selected. 

The following equations are used: 


Formula 21-1 


A—Availability 

U-Unavailability 

М-Ж of elements in subsystem 

k—# of elements required for system operation 
п-Ж of failures during time t (average) 
P(k,t)-Prob. of k failures during time t 
i-element, s-subsystem subscript 
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General formulas: 


Жы МТВЕ 
“ MTBF + MTTR 
U=1-A 


Series Subsystems: 





P(k,t) = Š 


n = t/MTBF 


=й (n) 
k! 








N > = 
А, = w Á, MTBF. = 1/5 (1/МТВЕ) MTTR. =} > ^: MTBF, 
1 1 s 
Parallel Subsystems: 
N A 
A, = 1 – T (1 = А9 МТВЕ, = I Ж MTTR, 
N 
MTTR, = 1/2 (1/MTTR,) 
1 
Parallel “К-о-М” 
MTTR, = MTTR/(N—-K-1) 
Е MTBE,| **" . (N—K)! (K-D)! 
MTBF, = MTBF, hé. Gm amm 
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INSTRUCTIONS 


Press the RUN key to start the program, then follow the prompt by entering a 
number between 0 and 6. 


0 — Start a new run 

1 — Enter MTBF and MTTR for a series subsystem 

2 — Enter MTBF and MTTR for a redundant parallel subsystem 

3 — Enter MTBF and MTTR for a redundant parallel k-of-n subsystem 
4 — Clear subsystem if a mistake was made during entry 

5 — [nclude a subsystem in an overall system 

6 — Find probability of x failures during time t 


EXAMPLE 


Problem Definition 


Consider a telephone office installation shown in Figure 21-1 consisting of a com- 
mercial and a stand-by power supply, three switching machines, two of which are 
required for operation, a group of incoming and a group of outgoing trunks con- 
nected to microwave antennas, all with known MTBFs and MTTRs. Find overall 
reliability parameters including probability of failure in 12 and 24 hours. 


Figure 21-1 
POWER SWITCHING INPUT/OUTPUT 
SUPPLY MACHINES (SERIAL) 
(PARALLEL ) (2-OF-3) 


COMMERCIAL 
POWER 


MTBF = 480 
MTTR=2 















INCOMING 
TRUNKS 


OUTGOING 
TRUNKS 








MTBF= 100 
MTTR=10 


STAND-BY 
POWER 


MTBF=5000 MTBF=4000 
MTTR=0.5 MTTR=0.5 





MTBF = IOO 
MTTR=2 
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PROGRAM LISTING 


OPERATION O - 67?2 


PARALLEL 

ENTER MTBF»MTTR?480:2 

A= 0.99585062 
МТВЕ= 480 

MTTR= 2 
OPERATION O - 672 

PARALLEL 

ENTER МТЕҒ»МТТК?100»10 

A= 0.99962278 
МТЕҒ- 4416.6667 
МТТК- 1.6666667 
OPERATION O - 675 

ASYS= 0.99962278 
МТЕҒ5Ү5- 4416.6667 
MTTRSYS- 1.6666677 


OPERATION O - 673 
K-OF -N 

ENTER МТЕҒ»МТТК?”100,2 
ENTER N»N?25»3 


A= 0.99880144 
МТЕҒ- 8535.33333 
МТТК- 1 
OPERATION O - 675 

ASYS= 0.99842467 
MTRFSYS- 701.9582 
МТТК5Ү5- 1.1061578 
ОРЕКАТТОМ O - 671 

SERTES 

ENTER МТЕҒ»МТТК?5000,0.5 
A= 0.99990001 
МТЕҒ- 3000 


МТТК- 0.50000025 
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OPERATION O - 671 


SERIES 

ENTER МТЕҒ»МТТК?4000»0.5 
A= 0.99977504 
A= 0.99977504 
MTRF- 2222.2222 
MTTR= 0, 50002885 
OPERATION O - 675 

ASYS= 0.97820007 
MTBFSYS= 332.93112 
МТТК5Ү5- 0.96097057 


OPERATION 0 - 67?6 
TIME:ND.FAILURES?24;,1 
F( 24y1 )=, 043050892 


OPERATION O - 676 
TIHEs;NO.FAILUREST12:5:1 
FC 12»12-,022015632 


OPERATION O - 676 
TIME:NO.FAILUREST24,2 
F(24522-,000269 376513 


Discussion of Results 


Looking at the subsystem first we see that the parallel power supply subsystem has 
a MTBF of 4416 hours, the 2-out-of-3 switching machine subsystem has a MTBF 
of 833 hours and the series incoming/ outgoing trunk subsystem has a MTBF of 
2222 hours. The overall system has a MTBF of 532 hours. To improve this factor 
would evidently require an improvement in reliability of the switching machines. 
The probability of one failure in 24 hours is 4.3%, of one failure in 12 hours is 
2.2%. The probability of 2 failures in 24 hours is only 0.097%. 


PROGRAMMING REMARKS 


The program fills nearly the whole 2K memory of the Timex/Sinclair 1000. For 
this reason the prompts had to be kept to a minimum. 
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PROGRAM LISTING 


15 GOSUE 800 

20 LET M=1 

40 LET N=0 

30 PRINT 

60 PRINT “OPERATION O - 6%"; 
70 INPUT Z 

80 PRINT Z 

90 GOTO 100Ж2410 
110 PRINT "SERIES" 
120 GOSUB 1000 
130 LET С=С+1/А 
140 LET D=D0KA/¢ ATE ) 
150 LET F=1/C 
160 LET H-Fx*C1-D)/TI 
170 LET Е=һ 
180 PRINT "A="sEy"MTBF="9Fs"MTTR=" ЭН 
190 GOTO 50 
210 PRINT "PARALLEL" 
215 GOSUR 1000 
220 LET С=С+1/В 
225 LET I-IXR/CAFEB) 
230 LET Е-1-П 
240 LET H=1/C 
290 LET F-EXH/TI 
260 GOTO 180 
510 PRINT "K-OF-N" 
515 GOSUE 1000 
520 PRINT "ENTER KsN7?"; 
525 INFUT K 
$30 PRINT К; 
335 INFUT К 
340 PRINT ">" R 
345 LET L-R-K 
450 GOSUB 900 
355 LET F=J 
360 LET L=K-1 
365 GOSUR 900 
370 LET F=FxXJ 
475 LET L-R 
380 С05ШЕ 900 
385 LET FzAXF/Jk*( А/К)ЖЖ‹ К-К) 


590 
595 
400 
410 
420 
430 
210 
320 
530 
540 
950 
560 
970 
610 
612 
614 
615 
620 
625 
630 
635 
640 
645 
650 
800 
810 
820 
700 
910 
920 
930 
940 
1000 
1010 
1020 
1030 
1040 
1050 
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LET Н-В/(К-К41) 
LET E-ZF/CF4H) 
GOTO 180 

GOSUB 800 

РКІМТ “SUBSYSTEM CLEARED" 
GOTO 50 

GOSUR 800 

LET M-MXE 

LET N=N+1/F 

LET Ғ-(1-М2/йМ/М 
LET X=1/N 


PRINT "ASYSz",.M»"MTRFSYSz" ,X»"MTTRSYSz",F 


GOTO 50 

PRINT "TIMEsNO.FAILUREST" » 
INFUT T 

РКІМТ Т; 

INFUT Q 

FRINT "+" yQ 

LET S=TXN 

LET L=Q 

GOSUR 900 

LET U=EXF -SkKSXKKQ/J 

PRINT "P("yTy" |" yQy" )=" sU 
GOTO 50 

LET C=0 

LET П=1 

RETURN 

(ЕТ J=1 

FOR 1=1 TOL 

LET J-JXI 

NEXT I 

RETURN 

PRINT "ENTER MTBF es MITRE" = 
INFUT A 

PRINT А; 

INPUT В 

PRINT "=" yË 

RETURN 
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MISCELLANEOUS 


22. Frequency of Musical 
Notes on the 
Well-Tempered Scale 


PROGRAM DESCRIPTION 


This program finds the frequency of a given musical note on the well-tempered 
scale. The scale assumes a constant frequency ratio of the 12th root of 2 between 
two adjacent notes. Johann Sebastian Bach was so excited by the idea of well- 
tempered scale that he wrote a series of piano compositions in the key of each note 
in that scale. Octaves in the program are numbered as negative or positive in- 
tegers, or 0, with note A in octave 0 having a frequency of 440-Hz. The frequency 
of notes doubles with each successively higher octave. 


INSTRUCTIONS 
Pressing the RUN key starts the program and issues two prompts for the note and 
for the octave number. The note should be entered as C, C+, D D+, Е, F, F+, 


G,G+,A, A+, orB; the octaveas — 2, —1,0,1,2...etc. As the Timex/ Sinclair 
1000 does not have a sharp (#) key, a (+) is substituted for a (#). 


EXAMPLE 


Problem Definition 


Find frequencies of A in octave 0, C# in the third octave, D in the — 4 octave, and 
C and B in the first octave. 
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Sample Run 


ENTER NOTE E.G. CCHI.. FÀ 
ENTER ОСТАУЕ?О 
A ОСТАЧЕ 0-440 HZ 


ENTER NOTE E.G. C;C D... PCH 
ENTER OCTAVE? 
Ct;OCTAVE 3=2217.4611 HZ 


ENTER NOTE E.G. C» Ct; Il, e PLI 
ENTER ОСТАУЕ?-4 
DrOCTAVE -4=18, 354048 HZ 


ENTER NOTE E.G. CrCtsIt.. TR 
ENTER ОСТАЧЕ? 3 
WRONG NOTE 


ENTER NOTE E.G. С>С%»П..?С 
ENTER OCTAVE?! 
CrsOCTAVE 1=523, 25113 HZ 


ENTER NOTE E.G. CsCtsll..?R 


ENTER OCTAVE? 
BsOCTAVE 1=987.7666 HZ 


Discussion of Results 
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The first example, as expected, yields 440-Hz, the fourth entry shows a wrong 
input (K) which yields an error return and repeat of the prompt. 


PROGRAMMING REMARKS 


Here, for a change, is a short program. Lines 100 and 135 perform the error check 
of inputs. If the note designation is not recognized as legitimate, then the variable 
I retains its assigned value of 100 and the ‘WRONG МОТЕ” remark is displayed. 


The frequency computation is performed in line 170. 
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PROGRAM LISTING 


10 
30 
S 2 
60 
70 
80 
70 
190 
110 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
1.52 
1.55 
150 
160 
170 
180 
190 


CLS 

PRINT 

PRINT "ENTER NOTE E.G. СЭС4»П,,?"; 
INPUT E$ 

PRINT В% 

PRINT "ENTER OCTAUE*?" ; 

INPUT B 

IF В<2ІМТ B THEN GOTO 80 

PRINT R 

LET 1=100 

IF R$-"C" THEN LET I=1 

IF R$-"Ct" THEN LET I-2 

IF R$-"I" THEN LET I=3 

IF R$-"Di" THEN LET I-4 

IF R$-"E" THEN LET I=5 
IF B$-"F" THEN LET I=6 
IF R$-"Ft" THEN LET I= 
IF R$-"G" THEN LET I=8 
IF B$-"Gt" THEN LET I=? 

IF B$-"A" THEN LET I=10 

IF R$-"At" THEN LET I-11 

IF B$-"B" THEN LET I=12 

IF І<>100 THEN COTO 170 

PRINT "WRONG NOTE" 

GOTO 50 

LET Ғ-2ЖЖ(( I-10 )/12+E )k44 0 

PRINT B$5",O0CTAUE "SBs"="s5F5s" HZ" 
GOTO 50 
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23. Artificial Intelligence 
Strategy Game 


PROGRAM DESCRIPTION 


This program, which is a take-off on the popular penny-matching HEADS/ 
TAILS game, is quite unique in that it learns by playing against the human oppo- 
nent. It thus demonstrates a strong intelligence trait. 

The game consists of the player calling either HEAD or TAIL. If the computer 
guesses correctly by making the same choice, then the computer wins. If the choice 
is different, then the player wins. If you, the player, play at random, then the com- 
puter also guesses at random and the number of wins and losses will be about the 
same. However, if you follow any of a number of simple strategies; e.g., play 
HEADS all the time, or play TAILS until you lose and then change to HEADS 
and vice versa, or play the same (HEADS or TAILS) if you lose, but alternate if 
you win, then the computer will outguess you after a few moves and beat you 
every time from then on, until you change your strategy. As soon as you change 
your strategy and the computer loses, it starts playing at random until it can lock 
on to your new strategy and win again. 

The program checks two previous moves in the game and consults its memory 
for the same type of situation played earlier. If a former situation is found, thena 
guess is made by the computer based on previous experience. If the player then 
plays as expected, the computer wins. If the situation did not occur earlier, the 
computer plays randomly and stores the outcome of its move with the player’s re- 
sponse to establish a new situation. 


INSTRUCTIONS 


Pressing the RUN key issues a prompt for a starter which initializes the RND 
BASIC function. The starter should be between 0 and 65535. Entering 0 will al- 
ways Start a different random sequence, otherwise the sequence is predetermined. 
The program then prompts for H(eads) or T(ails). It finally announces who won, 
shows the total score and continues with prompting for the next move. 


EXAMPLE 


Problem Definition 


The player (you) will first follow the strategy of alternating when losing, playing 
the same when winning. Then, when the computer catches up with your strategy, 
you will reverse your strategy. 
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STARTER?1 

(HOÓEAD OR CTOAIL'?H 

GAME 1sYOQU ШОМ; WIN COUNT=1 
(HOEAD OR € TOAIL?H 

GAME 2»YOU ШОМ; WIN COUNT=2 
(Н)ЕАП DR (TOAIL?H 

GAME 3»YOU LOST» WIN COUNT=1 
(HEAD OR CTOAIL?T 

GAME 4» YOU ШОМ», WIN COUNT=2 
НЕАП OR (TOAIL?T 

GAME 55.YOU WON» WIN COUNT=3 
(НЕАП OR CTOAIL?T 

GAME 65YOU LOST» WIN COUNT-2 
(НЕАП OR C TOAIL?H 

GAME 7,Ү00 LOST» WIN COUNT=1 
(НЕАП OR (ТОАП ?Т 

GAME 8»YOU LOST» WIN COUNT=0 
(Н)ЕАП OR (TOAILTH 

GAME 9,YOU LOST» WIN COUNT=-1 
(HOEAD OR ¢( TOAILTT 

GAME 10,Ү00 LOST» WIN COUNT--2 
(НЕАП OR (TOAIL?T 

GAME 11,Ү00 WON» WIN COUNT--1 
(НЕАП OR (C TOAIL?H 

GAME 12,Ү00 ШОМ, WIN COUNT-O 
(НЕАП OR (TOAIL?T 

GAME 135»YOU LOST» ШІМ COUNT--1 
(НЕАП OR CTOAIL?T 

GAME 14,Ү0/ LOST» WIN COUNT--2 
(НЕАП OR CTOAIL?T 

GAME 155YOU LOST» WIN COUNT--3 
(H)EAU OR CTOAIL?T 

САМЕ 16,/YOU LOST» WIN COUNT--4 
(HEAD OR ( TOAIL?T 

GAME 17,YOU LOST» WIN COUNT--5 
(НЕАП OR ( TIAIL?H 

САМЕ 18»YOU ШОМ» WIN COUNT--4 
(НЕАП OR (€ TOATL?H 

GAME 19,YOU LOST» WIN COUNT--5 
(НЕАП OR CTOAIL?T 

GAME 20,Ү00 LOST» WIN COUNT--6 
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(НЕАП OR € T)ñIL ?H 

САМЕ 21»YOU ШОМ, WIN COUNT--5 
(НЕАП OR € TDAIL'*H 

GAME 22,Ү00 LOST» WIN COUNT--6 
(HOEAD OR € TOAIL'T 

GAME 23, YOU LOST» WIN COUNT--7 
(НЕАП OR € TOAIL?H 

САМЕ 24»YOU LOST» WIN CDUNT--8 
(HOEAD OR C(CTOAIL? 


Discussion of Results 


You follow your first strategy through game 10. The computer is winning each 
game, starting with game six. When you change strategy at game 11 you are first 
winning but the computer catches up again starting with game 13. You then revert 
to your old strategy in game 18. After a few moves the computer catches up with 
you again. After 24 games you are down 8 points. 


PROGRAMMING REMARKS 


The outcome of various plays are stored in the dimensioned variable A. The check 
for correct input (H) or (T) is performed in line 200. The program requires the 
SLOW mode so that results can be continuously displayed during the INKEY$ 
sequence. 


PROGRAM LISTING 


10 CLS 

15 DIM ACS) 

18 SLOW 

20 PRINT “STARTER?” 5 

30 INPUT X 

40 PRINT X 

45 RAND (X) 

90 LET J=0 

60 LET I=0 

70 LET I=I+1 

80 IF I<3 THEN GOTO 140 
90 LET L-INT (C(MXNT1)0/2) 
100 LET W-4XGQt2XLtFt1 
110 IF ACU)-0 THEN GOTO 140 
120 LET 5-МЖА(С U) 
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130 
140 
170 
180 
185 
190 
200 
205 
210 
220 
250 
240 
250 
260 
270 
280 
290 
500 
510 
520 
330 
340 
550 
560 
570 
580 
590 
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GOTO 170 
LET Sz2XINT (2XRNI)0-1 


FRINT "C(HOEAI! OR CTOAIL?"; 
IF INKEY$-"" THEN GOTO 180 


LET T$-INKEY$ 


IF INNEY$-T$ THEN СОТО 190 
THEN GOTO 180 


IF T$z2"T" AND T$<>"H" 
PRINT T$ 
LET U=1 


IF T$-"T" THEN LET U--1 


LET Ј=Ј-5Ж0 

LET Zs52-I 

IF Z=0 THEN GOTO 320 
IF 7>0 THEN GOTO 340 


IF ACW)=0 THEN GOTO 310 


IF S=U THEN GOTO 320 
LET ACW )=0 

GOTO 320 

LET АЫ )= МЖО 

LET N-M 

LET Q=F 

LET M=U 

LET P=INT ((5Ж0412/2) 
LET АФ= "ЦОМ, " 


IF S=U THEN LET A$-"LOST;" 


FRINT "САМЕ "sIs",YOU 
GOTO 70 


" jas i" 


WIN COUNT="5J 


24. Multiple Dice Throw 


PROGRAM DESCRIPTION 


Many games; e.g., the popular YAHTZEE, and even computer modeling simula- 
tions, are based on the outcome of a multiple dice throw. This program permits a 
throw of several dice at once and displays the result of each multiple throw. The 
Timex/ Sinclair 1000 BASIC pseudo random generator RND in line 90, which 
does the ‘‘throwing,’’ works well in this application; however, if you have your 
own favorite, feel free to replace that line. 


INSTRUCTIONS 
Pressing the RUN key displays a prompt HOW MANY DICE for this throw. 
After a number n (one or more) is entered, the program displays n digits, where 


each digit is between one and six and corresponds to the spots on a die. Then the 
next prompt appears on the display for another throw. 


EXAMPLE 


Problem Definition 


Keep throwing dice, as shown in the listing. 


Sample Run 

HULTIPLE DICE HOW MANY DICE 72 
ёё 

HOW MANY DICE ?5 HOW MANY DICE 7?3 

4,2,2»3»5 29954 

HOW MANY ШІСЕ 73 HOW MANY DICE 77 

6,6,1 355»6,5»4,2,4 

HOW MANY DICE 71 HOW MANY DICE 75 

1 693962235 

HOW MANY DICE 71.5 HOW MANY DICE 76 

HOW MANY DICE 70 3?4»2»3»3,3 

HOW MANY DICE 710 

Irás 556,59 25156545 5 


99 


100 SCIENCE AND ENGINEERING PROGRAMS 


Discussion of Results 


The results are displayed. When we requested 0 or 1.5 dice, since both numbers are 
not acceptable, the prompt reappears. 


PROGRAMMING REMARKS 


Error check of the input is performed in line 70. The program is straightforward 
and short. Line 90 indicates how to transform the random number generated by 
RND from the range 0-1 into the range of integers between 1 and 6. 


PROGRAM LISTING 


10 CLS 

20 PRINT "MULTIPLE DICE" 

30 PRINT 

40 PRINT "HOW MANY DICE 7"; 
90 INPUT А 

60 PRINT А 

70 IF A<>INT A OR А<1 THEN GOTO 40 
80 FOR I=1 TO A 

90 PRINT INT (6XRNIO-T15; 

100 IF I=A THEN COTO 120 
110 PRINT ","; 

120 NEXT I 

150 GOTO 30 


25. Two-User 
Checking Account 


PROGRAM DESCRIPTION 


Here is a program to aid in keeping the financial side of a household running 
smoothly. In these days of two working spouses the problem frequently arises 
how to handle a joint checking account. Both parties may keep separate check- 
books in which they record separately all their transactions such as deposits, 
checks being written against the account, and transfers between the two spouses. 
When the bank sends a statement at the end of the month the following program 
will help in balancing the two checkbooks. The equations used are as follows: 


Spouse #2 balance = previous month’s balance of spouse #2 + 
deposits of spouse #2 — cleared checks of spouse #2 + 
transfers from spouse #1 to #2, also 


Spouse #2 balance = current bank balance — 
checkbook balance #1 — non-cleared checks #1 — 
earned interest (for NOW accounts) + bank charges. 


The spouse #2 balance is computed by the program in two ways and the results are 
compared. If they are the same then the account is balanced. If they are different, 
then there is an error either in the bank statement or, more probably, in one of the 
checkbook entries. 


INSTRUCTIONS 


Pressing the RUN key issues a prompt for one of the following commands: 

Z — Clears all registers, starts program 

X — Prompts for previous balance of spouse #2 

C — Prompts for each deposit to joint account by spouse #2 

V — Prompts for each cleared check by spouse #2 

B — Prompts for each internal transfer from account of spouse #1 to #2; use 
negative values for transfers from #2 to #1. 

N — Prompts for current checkbook balance of spouse #1 

M — Prompts for non-cleared checks written by #1 

A — Prompts for the current bank balance 

S — Prompts for earned interest (NOW checking accounts) 

K — Prompts for bank charges, if any 

D — Displays spouse 72 balance computed both ways 
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EXAMPLE 


Problem Definition 


The monthly bank statement just arrived showing a balance in the joint account of 
$1234.45. Spouse #1 has a checkbook balance of $751.23, and two checks written 
by #1 for $15.25 and $10.00 did not yet clear. The bank credited the account with 
$6.75 in interest. Spouse #2 during the month wrote three checks for $250.00, 
$125.00, and $50.00, but the last check did not yet clear. Spouse #2 had deposits 
credited for $150.00 and $175.00, and an internal transfer from spouse #1 for 
$300.00. The previous checkbook balance of spouse #2 was $203.22 and bank 
charges came to exactly $2.00. Find the current checkbook balance of spouse #2 
and see that the account is balanced. 


Sample Run 


2 FERSON CHECKING 


ENTER ZsXeCsVeBsNsMsAvSekK OR DPX 
PREV BAL NO, 27203,.22 


ENTER Z»X»C»U;E;N» Mss S» K OR LPC 
DEFOSITS NO. 271505; TOT=150 


ENTER Z»X»C»U:s Bs N» Mr 5 5» K OR DPC 
DEFOSITS NO. 27175; TOT=325 


ENTER Z»X»sC»U» EN» Ms As S» K OR DPP 
CL. CHECKS NO. 27250; TOT=250 


ENTER 2»Х»С»,У>Е»М»>М»Аэ5эК OR ПРИ 
CL. CHECKS NO. 27125; TOT=375 


ENTER ZeXeCsVeBReNeMsAsSsK OR ПРЕ 
XFERS 1 TO 27300,» TOT=300 


ENTER ZsXsC US ENM As SK OR ТРМ 
BAL М0.17751.23 


ENTER ZsXsCoVseBsNotsAsSek OR ОРМ 
OUTST CHKS М0.1715.25», ТОТ-15.25 


ENTER 2»Х>С»,У»,В»М»М»>йһ»,5эК OR ПРМ 
OUTST CHKS М0.1710>,Т0Т-25, 25 
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ENTER 2»Х>С>У,Е ЭМ, M ñ | S5» K OR DPA 
BANK R6L71254.45 


ENTER 7»,Х,С»У, КУМ Mrs 5» K OR LPS 
INTEREST?T6.73 


ENTER 7»Х:С,У»,Е»,М»М»,А»б,К OR [PK 
BANK CHARGES?T2:s TOT=2 


ENTER ZsXeCeVeBReoNestsAsSeK OR ПРП 
BAL NO. 2=453,22,453,22 
ACCOUNT BALANCED 


Discussion of Results 


The account fortunately balanced; both computations led to the same result. 
Notice that the $50 check written by spouse #2, which did not clear yet, did not 
have to be entered. 


PROGRAMMING REMARKS 


Branching is performed via the INKEYS function and input is then compared with 
the recognized commands Z, X, C... Line 580 takes care of the roundoff error of 
the Timex/Sinclair 1000. Notice that in the subroutine starting in line 1000 all 
variables are set to 0. Why not use the CLEAR command instead? The reason is 
that then statements such as those in line 190 or 250 would return error 2 (unde- 
fined variable). 


PROGRAM LISTING 


10 CLS 

15 РКІМТ "2 PERSON CHECKING" 
20 SLOW 

30 GOSUR 1000 

40 PRINT 

45 PRINT "ENTER ZsXeCeVsBsNeMsAsSeK OR LP's 
50 IF INKEY$-"" THEN GOTO 50 
3o PAUSE 5 

60 LET A$-INKEYS$ 

65 IF А%-ІМКЕҮ% THEN GOTO 65 
70 PRINT f$ 

90 IF А%<>"7" THEN GOTO 110 
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95 
100 
105 
110 
120 
150 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
500 
510 
520 
330 
340 
550 
360 
370 
380 
390 
490 
410 
420 
430 
440 
450 
455 
460 
470 
480 
490 
200 
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GOSUR 1000 

PRINT "SYSTEM СІҒАКЕП" 
GOTO 40 

IF А%<>"Х" THEN GOTO 160 


FRINT “PREV BAL NO. 27?"; 
INFUT X 
PRINT X 
GOTO 40 
IF А%<>"С" THEN GOTO 220 


PRINT "DEFOSITS NO. 27"; 


INPUT Р 

LET С=С+Р 

PRINT FPs"#TOT=" +C 

GOTO 40 

IF A$<>"V" THEN GOTO 280 
PRINT "CL. CHECKS NO. 27" 
INFUT Y 

LET U=U+Y 

PRINT Ys" »TOT-"sV 

GOTO 40 

IF А%“>"Е" THEN GOTO 340 
FRINT "XFERS 1 TO 2?"; 


INFUT L 

LET В-ЕН 

PRINT Ls" ,TOT-"sE 

GOTO 40 

IF А%<>"М" THEN GOTO 390 
РКІМТ "BAL NO0.1?"; 

INFUT N 

PRINT N 

GOTO 40 

IF А%<>"М" THEN СОТО 450 
PRINT “OUTST CHNS NO0.17?"; 
INPUT J 

LET MzM4J 

PRINT Js" » TOT" yM 

GOTO 40 

IF 6$2-"A" THEN GOTO 490 
PRINT "BANK ВАГТ", 

INFUT А 

PRINT А 

GOTO 40 

IF А%<>"5" THEN GOTO 540 
PRINT "INTEREST?" ; 


510 
520 
930 
340 
550 
260 
370 
380 


290 
610 
620 
630 
640 
650 
655 
660 
670 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
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INFUT S 

FRINT S 

GOTO 40 

IF А%<>"П" THEN GOTO 610 

LET F=X+C-V+tk 

LET G=A-N-M-S+tK 

PRINT "BAL NO. 2="yFy" |" yG 
IF ABS (F-G)<.001 THEN PRINT 
"ACCOUNT BALANCED" 


GOTO 40 

IF А%<>"К" THEN GOTO 660 
PRINT "BANK CHARGES?" 5 
INFUT F 

LET K=K+F 

PRINT Ғ%“»ТОТ-"%К 

GOTO 40 

PRINT "REPEAT" 

GOTO 40 

LET C-0 

LET V=0 

LET В=0 

LET M=0 

LET К=0 

LET Х=0 

LET А=0 

LET М=0 

LET S=0 

RETURN 
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APPENDIX 


Programming Hints 
and Techniques 


While writing and debugging programs in this book I came up with a number of 
programming hints and techniques that I would like to share with other Timex/ 
Sinclair 1000 users to save them some grief and frustration. They are listed here. 


1. Roundoff 
If you deal with floating point numbers do not make a decision on equalities of 
twonumbers(e.g., IF X = YTHEN.. .). The TS 1000 handles roundoff strangely; e.g., 


100 LET X =3.25 —3 – .25 
110 PRINT X 


will display 1.164 ... Е- 10. Therefore always take the absolute value of the dif- 
ference and compare it with a small number. 


2. Running Out of Memory 

When you write a long program which causes the Timex/Sinclair 1000 to run 
out of memory, the computer will start dying ‘‘gracefully’’ and will make your life 
miserable. One of the first symptoms of memory ‘‘tightness’’ is when the com- 
puter ignores the EDIT function (Shift 1). However, as the Timex/Sinclair 1000 
shares display and program memory you can get around this problem. For example, 
if you want to edit line 350 while you are running out of memory, do the following: 


LIST 350 Listing appears on the screen with line 350 on top 
CLS Screen goes blank 
EDIT Presto, Line 350 appears on the bottom of the screen ready to be edited 
3. Multiple Branches 

The 1000 lacks a computed GOTO. To get around this deficiency, number 
various branches of your program in multiples of 100 and use the statement 
GOTO 100*A 


The variable A should of course be an integer. 
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4. Leaving a Loop 

If you are stuck during program execution in a tight loop with an INPUT state- 
ment in it, you may not be able to get out of the loop without unplugging the com- 
puter. The BREAK key will simply be interpreted by INPUT as a space. If the 
expected input is numerical you can still get out of the loop by keying in, for exam- 
ple, Shift Enter (Function) LN 0, and getting an ‘‘A”’ error return. If the expected 
input is alphanumeric, then your only hope is to try to fill the screen and to get out 
with a <5”” error return. 


5. INPUT Statement 

Precede all INPUTS with a prompt. An INPUT statement does not produce a 
question mark automatically as in some other BASICSs. Follow an INPUT state- 
ment with a PRINT to display what you have just entered. Otherwise you will 
never know. A typical sequence could read as follows: 


10 PRINT “ENTER Х?”; 
20 INPUT X 
30 PRINT X 


6. Use of INKEY$ 

INKEYS3 is used as a convenient means for program branching. Within a loop 
with INKEY$ the computer should be in the SLOW mode or the program should 
have PAUSE statements, otherwise the screen will remain blank with no visible 
prompts. Here is a foolproof sequence to detect, for example, the **X"' key. It will 
operate in both SLOW and FAST mode: 


100 PRINT “УУНАТ NEXT?"'; 

110 PAUSE 10 

120 IF INKEY$= '' ” THEN GOTO 110 
130 LET A$ = INKEY$ 

140 PRINT A$ 

150 PAUSE 10 

160 IF A$ = INKEY$ THEN GOTO 150 
170 IF A$ < >“Х” THEN GO TO 100 
180 REM “X” HAS BEEN DETECTED 


7. Use of Quotes 
To indicate an alphanumeric null use Shift P (**) twice and not Shift Q (** ’’). 
Shift Q can only be used to display quotes inside alphanumeric strings. 


8. Variable Initialization 

Every variable has to be initialized with an assignment or a DIM statement be- 
fore it can be put on the right side of an assignment statement.Otherwise the pro- 
gram execution will stop with a “2”? error. Though annoying at first, this feature is 
great for program debugging. It will flag any misspelled or uninitialized variables. 
But beware of the CLEAR command which will obliterate all variable initializa- 
tions (assignment or DIM). In fact I have found very little use for the CLEAR 
command. 
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9. Dimensioning of Strings 

Be careful when using dimensioned strings. A singly dimensioned string can be 
only one character long. The second dimension is interpreted by the 
Timex/Sinclair BASIC as the length of the string. 


10. Generating Similar Program Statements 

A fast way to generate several similar statements in your program is to write the 
first statement, then edit it and also edit the program line number. Editing a pro- 
gram line number generates two statements — one with the old program line 
number, one with the new one. The Timex/Sinclair’s BASIC permits editing of 
line numbers, while most other BASICs do not. 


11. Deleting the Input Buffer 

When entering data or programs from the keyboard the characters are tempo- 
rarily stored in a buffer and are displayed on the bottom of the screen. Pressing 
the ENTER key transfers the buffer into main memory. If at some point you de- 
cide to erase the current buffer you may do it quickly by pressing Shift 1 (EDIT) 
rather than deleting it character by character or token by token with the Shift 0 
(DELETE) key. 


spectrum Conversions 


The ZX Spectrum is a further step up the ladder of the Sinclair’s home computer 
family. It retains all the good features of the ZX-81 and the Timex/ Sinclair 1000 
and adds many new ones. The ZX Spectrum computer in its basic configuration 
comes with 16K of RAM instead of 2K in the Timex/Sinclair 1000, has color 
and extended graphic capabilities, and adds several new BASIC functions. For- 
tunately most ZX-81 and Timex/Sinclair 1000 BASIC functions are identical to 
those featured in the ZX Spectrum. Therefore, most of the programs featured in 
this book will run on the ZX Spectrum exactly as shown, without any modifica- 
tions. 

However, for better readability all variable names shown in this book in upper 
case should be entered in the ZX Spectrum in lower case, its standard operating 
mode. Certain functions have been slightly changed in the ZX Spectrum in func- 
tion or in name. For example, the RAND function used in program #14 and #23 is 
called RANDOMIZE. The SLOW and FAST functions used in several pro- 
grams in this book do not exist on the ZX Spectrum and should be omitted. 
Another minor change is replacement of the exponentiation sign (**) by an up- 
arrow. One difference between the two computers which should be carefully 
watched is the change in values of the CODE function. The CODE function in the 
ZX Spectrum is very similar to the ASC function found on most other BASIC 
computers. For proper operation the following changes should be made in pro- 
grams #6 and #8 which make use of the CODE function. 


Program #6. Hex-to-Decimal and Decimal-to-Hex Conversion. 


Enter hexadecimal digits A-F in upper case. Replace lines 260 and 510 and add 
lines 261, 262, and 511 as follows: 


260 LET с = CODE а$ (1-1) 

261 IF c < = 57 THEN LET c = c — 48 

262 IF c > 57 THEN LET c = c — 55 

510 IF e < = 9 THEN LET b$(k) = CHRS$(e + 48) 
511 IF e > 9 THEN LET b$(k) = CHR$(e + 55) 


Program £8. Relative Offset Computations for Machine Code Programming. 


Enter hexadecimal digits A-F in upper case. Replace lines 280, 290, and 370 and 
add lines 281, 291, and 371 as follows: 


280 IF g <= 9 THEN LET e$ = CHR3(g + 48) 

281 IF g > 9 THEN LET e$ = CHR + (в + 55) 

290 IF g < 9 THEN LET w$ = СНК$(2 — 16 * g + 48) 

291 IF g > 9 THEN LET w$ = CHR3$(z — 16 * g + 55) 

370 IF CODE c$(l-i) < = 57 THEN LET г = CODE c$(l-i) — 48 
371 IF CODE с%(1-і) > 57 THEN LET г = CODE c$(1-i) - 55 
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arrival rate, 80 
assembler, code, 31 
attenuation, 7 
availability, 85 


bank balance, 101 

base conversion, 22, 25 
binary, base, 22 
binomial, coefficients, 44 
bit error rate, 18 

buffer, deleting, 108 


capacitance, 1 
checking account, 101 
coil, parameters, 1 
combinations, 44 
complex, numbers, 72 
cryptography, 28 
curve, fitting, 35 
curve, French, 35 


data, link, 18 

data, transmission, 18 

dice, throw, 99 

dielectric, constant, 12 
distribution, Gaussian, 48, 56 
distribution, binomial, 51 
distribution, exponential, 80 
distribution, hypergeometric, 51 
distribution, normal, 48, 56 
distribution, uniform, 56 


equations, differential, 65 
equations, transcendental, 69 
Eratosthenes, sieve, 28 

error, function, 18, 48 
extrapolation, 35 


facility utilization, 80 
factorials, 44 

failure rate, 85 

Fourier, coefficients, 16 
fractions, binary, 77 
fractions, decimal, 77 
function, even, 16 
function, odd, 16 
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games of chance, 99 
goto, computed, 106 


hexadecimal conversion, 22, 25 
hints, programming, 106 


impedance, characteristic, 12 
impedance, matching, 12 
inductance, coil, 1 

initial conditions, 65 
initialization, 107 

inkey$, 107 

intelligence, artificial, 95 
interpolation, 35 


Lagrange polynomial, 35 


machine language, 31 


mean time before failure (MTBF), 85 


mean time till repair (MTTR), 85 
mean, distribution, 56 
memory, ''tightness"', 106 
microprocessor, 31 

modulation diff, phase shift, 19 
modulation, OOK, 18 
modulation, amplitude shift, 18 
modulation, frequency shift, 18 
modulation, phase shift, 19 
modulation, pulse shift, 18 
Monte-Carlo simulation, 56 
musical tones, 92 


newton, iteration, 69 
number theory, 22 


octaves, 92 


pad, matching, 7 

pad, minimum loss, 7 
penny matching, game, 95 
periodic function, 60 
permutations, 44 

plotting, 35, 40 
polynomial, 35 

prime factors, 28 

pseudo random deviates, 56 
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queuing equations, 80 
quotes within quotes, 107 


redundancy, 85 
references, 109 

registers, arithmetic, 72 
relative addressing, 31 
relative offset, 31 
reliability, 85 

response time, 80 
reverse Polish (RPN), 72 
roundoff, 106 
Runge-Kutta, 65 


sampling, universe, 51 
scaling, readable, 40 
servers, 80 

service rate, 80 

shielding, 12 
signal-to-noise, 18 
smoothing, function, 35 
stack (RPN), 72 
standard deviation, 48 
statement number, editing, 108 
strategies, play, 95 
string, dimensioning, 108 


tempered scale, 92 
transmission, data, 18 
transmission, line, 12 
utilization, 80 


waiting time, 80 


YAH-T-ZEE, 99 


General Interest 
Computers 





Twenty-five science and engineering programs for the 
Timex/Sinclair 1000 bring professional sophistication 
їо the world of low-cost computing! 


With detailed documentation and liberally illustrated 
by sample runs, these programs solve the most fre- 
quently encountered problems in the areas of: 


electrical engineering 
data transmission 

` number theory 
computer programming 
computer-generated plotting 
probability and statistics 
mathematics | 
operations research 





These compact, clearly-written programs are compati- 
ble with the Timex/Sinclair Spectrum computer, and 
can easily be translated to run on other BASIC 
computers. 
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